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ABSTRACT 


A  series  of  experiments  were  conducted  on  the  influence  of  an  electric  field  on  free 
convection  heat  transfer.  A  heated  plate  was  mounted  in  a  Mach-Zehnder  interferometer, 
and  the  effects  of  applied  electrical  fields  were  studied.  Large  changes  in  effective  heat 
transfer  were  found.  Analytical  studies  were  made  of  this  phenomenon.  A  suitable  mathe¬ 
matical  model  was  formulated  on  the  basis  of  the  Von  Karman  integral  equations  and  the 
electrostatic  field  equations.  Correlation  between  the  interferometer  test  data  and  the 
mathematical  model  was  achieved. 

Exploratoiry  tests  of  fluid  orientation  and  motion  under  the  influence  of  dielectric  forces 
were  also  conducted,  -an^he  potential  usefulness  of  the  phenomenon  is  graphically,  por¬ 
trayed  in  pictures  showing  the  motion  of  water  droplets  in  a  “Zero-G”  model  when  a  field 
is  applied. 
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•  .  ^  •  Po 

ambient  pressure 

Ib/sq.  ft.  . 

•  :  Pc. 

centerline  stagnation  pressure 

Ib/sq.  ft.  '•  .  •  •  ■ 

■  ••  '  % 

internal  heat  generated 

BTU/hr-cii.  ft. 

•  **r- 

energy  radiated 

BTU/hr-cu.  ft.  •  '  . 

•  •  R'-- 

radius  of  electrode 

feet  ;  ■  ■  ■ 

••  .  '  /■ 

distance 

meters  i 

xli 

*  ^  '  ’  *  •  *  ^  •  j 

. 
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LIST  OF  SYMBOLS  (Cortt'd) 


SYMBOL 

DEFINITION 

UNITS 

Ra 

distance 

meters 

r 

radial  distance 

meters 

radius  at  edge  of  intense  corona 

meters 

•  T 

temperature 

•f..  • 

o 

ambient  temperature 

'F  . 

T 

w 

wall  temperature 

•f  . 

t 

transformation  plane 

dimensionless 

U 

impinging  velocity 

ft./8ec 

u 

velocity 

ft./sec 

u 

X 

velocity  in  boundary  layer 

ft./sec  ■ 

U' 

(» 

velocity  outside  of  boundary  layer 

ft./sec 

V  • 

applied  potential 

volts 

voltage  at  start  of  corona  current 

volts 

V 

.  velocity  '  •  ' 

■  ft./sec  ' 

total  ion  velocity  ’  ' . 

ft./sec 

w 

complex  number 

X 

body  force 

Ibs/cu.  ft.  ■ 

X 

coordinate  distance -distance  along  plate 

inches 

•  y 

coordinate  distance-distance  perpendicnlar 
to  plate  .  .  • 

•  inches 

z 

z  s  x/a  dimensionless  distance  along  pla^. 

dimdnSioiiiess 

a 

complex  number  s.x  +  ly  • 
thermal  diffusivity  k/pCp  '  ’ 

ft.Vhr' 

)3; 

stagnation  pressure  constant  ■  ' 

1/sec 

y 

pressure  ratio  parameter 

•  dirndnsionless' 

L 

5U  +  ^ 

3x  ay  9Z  ■ 

s’ 

velocity  boundary  layer  thickness 

feet  ■ 

6o 

boundary  layer  thickness  at  origin 

.  feet  • 

6t 

thermal  boundary  layer  thickness  ’  •  . 

feet 

ratio  of  thermal  to  velocity  bountopy  layer 
thickness 

dimensionless 

€ 

permittivity . 

farads/meter 

LIST  OF  SYMBOLS  (Cont’d) 
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SYMBOL 

C 

V 
$ 

Be 

X 

ti 

ua 

V 

P 

«: 

(Tc 

T 

« 

\li 


DEFINITION  UNITS 


transformation  plane 

heat  transfer  variable,  defined  in  text 
temperature  difference,  T-T^ 

temperature  difference, 

T 

angle  subtended  by  electrode 

radians 

charge  per  unit  length 

coulombs/meter 

dynamic  viscosity 

lb/ft.-hr 

current 

microamperes 

kinematic  viscosity 

ft.® /hr 

density 

lb  /cu.  ft. 
m 

charge  density 

coulombs/cu.  meter 

electrical  ccmductivity 

1/ohms 

time 

seconds 

dissipation  function;  defined  in  text 

stream  function,  defined  in  text 
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INTRODUCTION 

The  influence  exerted  by  an  electric  field  on  gases  and  liquids  has  been  known  for  many 
years.  Primary  attention  has  been  given  to  studying  electrolytic  phenomena  and  gas  dis¬ 
charges.  however,  and  relatively  little  work  has  been  done  on  the  effects  of  electric  fields 
acting  alone  on  the  flow  of  fluids.  Such  fleld  effects  tend  to  be  rather  subtle  in  contrast  to 
the  relatively  large  fluid  body  force  generated  by  the  interaction  of  an  electric  field  and' 
a  transverse  magnetic  field.  This  latter  fluid-electromagnetic  interaction  gives  rise  to 
the  whole  field  of  magnetohydrodynamics  (MHD).  The  interest  in  'MHD  has  received  wide 
attention  in  recent  years  because  of  potential  applications  in  power  generation,  controlled 
fusion,  and  space  propulsion.  The  subtle  electric  field  effects,  however,  merit  considera¬ 
tion  aj)art  from  the  mass  of  the  work  done  on  magnetohydrodynamic?.  The  fluid-electric 
interactions  are  generally  distinct  from  those  in  the  MHD  field. 

An  electric  field  can  influence  not  only  fluids  containing  charged  particles  arid  conduct¬ 
ing  fluids,  but  neutral  nonconductii^  fluids  as  well.  This  wide  range  of  action  may  provide 
controllable  body  forces  within  the  fluid  without  the  high  degree  Of  iom2:ation  customarily 
used  in  magnetohydrodynamics.  Changes  in  fluid  properties,  such  as  viscosity  or  con-.  ■ 
ductivity,  can  also  occur  in  the  presence  of  a  field,  and  thus  the  entire  behavior  of  fluid 
flow  and,  heat  transfer  may  be  affected. 

It  is  the  purpose  of  this  research  to  explore  the  possible  interactions  tetween  electric- 
fields  and  fluids.  Because  of  the  wide  range  of  possible  influences,  an  extensive  survey 
of  possible  electric  field-fluid  interactions  was  conducted  to  provide  a  sound  basis  for 
the  specific  areas  of  research  chosen.  A  summary  of  the  interaction?  presented  in  Ref-  . 
erence  1,  includes  the  basic  electrical  phenomena  involved,  descriptions ‘of  the  various 
interactions  observed  in  experiments,  and  a  discussion  of  the  potential  application  of  the 
phenomena.  .  .  •  ^ 

Two  specific  electric-fluid  interactions  were  Selected  for  these  more  detailed  investi¬ 
gations:  the  effects  of  ion  motion  in  gases,  and  the  motion  of  neutral- particles.  In  the- 
study  of  the  ion  motions,  a  combined  experimental  and  analytical  program  of  the  effects 
of  a  corona  discharge  on  the  free  convection  from  a  heated  plate  was  undertaken.  A  Mach-  • 
Zehnder  interferometer  was  used  to  record  results  of  experiments,  and  the  analytical 
approach  was  based  upon  the  Navier  Stokes,  energy.  Von  Karman,  and  electric  field 
equations.  Correlations  between  the  theoretical  solutions  and  the  test  ’data  were  accom¬ 
plished.  The  motion  of  neutral  particles  as  represented  by  dielectrophoresis  was  studied 
only  in  an  exploratory  fashion  to  verify  the  possible  actions  which  could  be  obtained  . 
through  the  application  of  an  electric  field.  An  explanation  of  the  corona'  wind  and  cUelec- 
trophoresis  is  included  in  Reference  1. 

EXPERIMENTAL  INVESTIGATION  OF  THE  EFFECT 
OF  CORONA  WIND  ON  FREE  CONVECTION 

This  section  covers  the  experimental  work  done  in  studying  one  possible  interaction.- 
between  an  electrical  field  and  a  fluid.  The  experiment  chosen  is  that  of  a  corona  dis¬ 
charge  impinging  on  a  heated  flat  ^ate.  The  section  will  describe  the  experimental- 
phases  of  the  work  and  a  description  of  the  phenomena  observed.  • 


Manuscript  released  by  the  author  on  June  1962  for  publication  as  an  ASD  Technical 
Documentary  Report. 


1 


o 


ASD-TDR-62-650 

AIMS  OF  THE  CORONA  WIND-HEAT  TRANSFER  TESTS 

In  order  to  study  the  possible  effects  which  could  be  achieved  by  electric  fields  actii^ 
on  ions  in  gases,  it  was  necessary  to  chose  a  phenomenon  which  could  provide  a  suitable 
source  of  ions.  Many  possible  sources  exist  including  combustion  flames,  ultra-violet 
radiation  of  the  gas,  and  others.  A  controllable  source,  and  one  which  is  simple  to  obtain, 
is  that  of  the  ions  in  a  corona  discharge  (Ref.  2). 

The  phenomenon  of  corona  wind  is  discussed  in  detail  in  Reference  1,  but  the  basic 
action  will  be  reviewed  briefly  here.  If  a  fine  needle-point  electrode  is  located  a  short 
distance  away  from  another  blunt  electrode,  a  stream  of  air  moves  away  from  the  point 
when  high  voltages  are  applied  between  the  point  and  the  other  electrode.  This  is  illus« 
trated  in  Figure  1. 


Figure  1.  Corona  wind  from  a  charged  point 


Around  the  fine  corona  point  an  intense  electrical  field  exists.  Ions,  existing  in  air  natuzaMy 
or  created  by  this  intense  field,  are  accelerated  by  the  field.  These  ions  collide  with  neu¬ 
tral  molecules  of  air  and,  near  the  point,  further  ionization  takes  place.  Away  from  the 
region  of  the  point,  the  collisions  are  elastic  and  a  streaming  of  the  air  results.  The  net 
effect  is  the  corona  wind  moving  away  frcxn  the  charged  point.  The  velocities  obtainable 
from  the  corona  wind  in  air  are  very  small,"  of  the  order  of  a  few  feet  per  second.  If  this 
slight  streaming  were  made  .to  impinge  upon  a  heated  plate,  however,  a  change  In  the  flow 
field  and  temperature  distribution  should  result.  The  free  conVectlCMi  from  the  heated  plate 
in  still  air  would  also  provide  a  very  sensitive  phenomenon  to  exjjlore  with  the  corona 
wind  or  with  any  other  electrostatic  influence  which  nilghi  exij^.  Cc^ise^uentlyi  a  com¬ 
bined  free  conyection-eiectrostatic  test  was  planned. 

To  actually  conduct  such  a-  test,  a  Mach  Zehnder  Interferometer  was  utilized  (Ref.  3). 
Using  the  Interference  patterns  of  light,  the  interferometer  can  record  very  small-changes  ;• 
In  the  temperature  field  around  the  test  specimen.  Any  changes  Induted  by  the  corona  wind 
or  electrostatic  field  should  thus  be,  recorded  by  the  interferometer  as  changes  in  the  inter¬ 
ference  pattern.  Thus  the  experiment  actually  conducted  is  based  upon  studying  a  very  sen¬ 
sitive  phenomenon  (free -convection)  with  an  extremely  sensitive  instrument  (the  interfero¬ 
meter)  to  detect  any  potentially  small  influences  of  the  applied  electrostatic  field.  The  data 
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obtained  from  the  interferometer,  although  erf  considerable  value,  pertain  only  to  tern-  . 
perature  distributions.  Since  no  Information  could  be  obtained  with  this  instrument  about 
the  velocity  distributiem  in  the  vicinity  of  the  plate,  pressure  effects  were  studied  separa¬ 
tely  and  the  velocities  were  then  deduced. 

Many  electrode  configurations  were  conceived  to  provide  a  suitable  corona  discharge 
and  the  field  shaping  necessary  to  control  the  phenomena.  Electrodes  that  would  provide 
a  field  without  a  corona  discharge  were  constructed  and  used  to  provide  base  line  data 
with  no  ions  and  to  explore  any  possible  effects  not  directly  associated  with  ion  motion. 
Because  the  corona  wind  tests  indicated  some  very  interesting  results,  most  of  the  work 
was  concentrated  ufxxi  the  corona  discharge.  The  work  with  the  other  electrode  configura¬ 
tions  were  exploratory  only. 

DESCRIPTION  OF  EXPERIMENTAL  APPARATUS 

Heat  Transfer  Equipment 

A  Mach-Zehnder  interferometer  was  used  to  study  the  effects  on  heat  transfer.  The 
theory  of  operation  of  an  interferometer  is  covered  in  numerous  publications  and  will  not 
be  discussed  in  detail  (Refs.  3,  4).  Both  infinite  fringe  and  finite  fringe  patterns  were  used 
in  the  study.  Because  the' infinite  fringe  patterns  gave  graphic  and  readily  usable  results, 
most  of  the  runs  were  made  with  this  setting.  The  phenomena  were  observed  both  with 
monochromatic  light  of  5641  A  and  with  white  light.  All  data  were  taken  with  the  mono¬ 
chromatic  light  source,  although  the  white  light  was  effective  in  visually  studying  the 
phenomena.  '  .  .  • 

Figure  2  shows  the  over-all  test  apparatus  with  the  interferometer  loqated  in  the  back¬ 
ground.  Depicted  are  the  high  voltage  supply  at  the  left,  the  working  sectirai  of  the  inter-  ' 
ferometer  in  the  center,  and  the  free-cohvection  heater  power  supply  on  the  table.  Figure 
3  indicates  the  Polaroid  Camera  located  at  the  rear  of  the  interferometer,;which.was  ' 
used  to  record  the  fringe  patterns  obtained.  The  controls  at  the  right  are  used  to  adjust' 
the  interferometer  to  the  fringe  pattern  desired. 

The  interferometer- was  located  in  a  small  building  isolated  from  outside  traffic  and  . 
disturbances.  All  windows  were  covered  oyer  and -taped  to  keep  out  extraneous  light.  ■ 
Since  random  convective  currents  in  the  room. caused  occasional  trouble  In. using  the , 
interferometer',  all  heat  was  turned  off  in  the  room  for  several  hbUrs  pirior  td  testing  eiid 
throughout  the  tests.  .Likewise,  even  though  the  room  had- been  apparently’sealed,  changes- 
in  atmospheric  conditions,  such  as  strong  wind,  bright  sunlight,  or  rapidly  changing  tem¬ 
perature,  caused  fluauations  in  the  action  of  the  interferometer.  Consequently  testing 
was  accomplised  only  on  relatively  calm  days  with  all  heat  turned  oft  ' 

The  specimen  used  to  study  free  convection  was  an  electrically  heated  flat  aljimlnum  - 
plate  6  inches  by  10  inches  In  size.  In  heat  transfer  tests  conducted.preyiously  by  per-  ' 
Bonnel  of  the  Air  Force  Institute  of  Technology,  the  temperature  distribution  of  the  plate  • 
had  been  found  to  be  very  uniform  over  a  range  of  temperatures.  The  plate  was  heated 
by  electrical  resistance  heatets  inbedded  in  the  rear  of  the. plate  and  well  Insulated  elec¬ 
trically  from  the  metal  place  itself,  Iron-constantan  thermocouples  were  located  In  the 
plate  in  the  positions  shown  on  Figure  4.  The  thermocouple  junctions  were  located  in  holes 
.drilled  in  the  plate  and  later  filled  in,  so  that  the  front  of  the  . plate  was.  smooth.  During 
the  course  of  the  tests,  four  reference  temperatures  on  the  plate  were  read  whenever  data 
were  taken  with  the  Interferometer.  These  temperatures,  remained  within  l*F  of  each  other. 
,  except,  in  a  few  cases  in  which  the  difference  was  slightly  over  1*F, 
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The  plate  is  shown  mounted  venlcally  in  the  interferometer  in  Figure  5.  Shown  also  in 
the  figure  are  a  typical  electrode  and  the  traversing  mechanisms  used. 

Temperatures  were  read  by  means  of  a  Rubicon  Instruments  potentiometer  using  an  ice- 
water  tath  reference  temperature  for  the  thermocouples. 

A.C.  current  was  used  to  power  the  heaters  in  the  plate  in  the  early  phases  of  the  test. 
The  power  dissipated  in  the  heaters  was  controlled  by  means  of  voltage  regulators.  The 
variation  of  temperature  over  the  plate  could  be  controlled  by  adjusting  the  applied  volt¬ 
age  to  specific  heaters,  and  the  mean  plate  temperature  could  be  varied  by  a  master  volt¬ 
age  control.  Although  this,  power  supply  provided  very  good  temperature  regulation  for 
the  free  convection  tests  without  applied  fields,  it  was  found  to  cause  difficulty  with  the 
microammeter  used  to  measure  corona  currents.  During  initial  studies  of  the  combined 
corona-heat  transfer  tests,  the  microammeter  gave  very  erratic  readings  whenever 
the  A.C.  heater  supply  was  turned  on.  Because  of  this,  the  A.C.  power  supply  was  re¬ 
placed  by  automotive  storage  batteries  to  provide  a  ripple-free  D.C.  power  source. 
Variations  in  voltages  required  for  temperature  control  were  secured  by  a  suitable  rheo¬ 
stat.  This  arrangement  eliminated  the  difficulty  and  the  microammeter  was  steady  through 
all  its  ranges.  Plate  temperatures  were  adjusted  to  uniform  values  and  remained  within 
the  1*F.  tolerance.  Figure  6  shows  the  heater  power  supply.  The  D.C.  supply  controls  are 
at  the  front  of  the  table  and  the  A.C.  Varlacs  are  at  the  rear  of  the  table.  TTie  batteries  are 
below  the  table,  and  the  voltmeters  and  ammeters  to  monitor  heater  input  power  are  in 
the  center  of  the  table.  The  thermocouple  potentiometer  is  in  the  rear  of  the  table  and  a 
mlcroamnieter  used  to  measure  corona  current  is  at  the  right, : 

Electrical  Apparatus 

A  0'-50,000  volt  power  supply  was  used  to  provide  the  electrical  field  for  the  corona  dis¬ 
charge.  It  had  a  current  capacity  of  up  to  2  milllamperes.  which  was  more  than  adequate 
for  the  tests  accomplished.  The  design  of  the  power  supply  allowed  .only' a  negative  ground; 
because  of  the  arrangement  of  the  test  equipment  only  a  positive  point  corona  could  be 
studied.  A  Grelbach  microammeter  used  to  measure  current,  had  ranges  from  0-1  micro¬ 
amperes  up  to  0-1000  microamperes. 

Three  instruments  vvere  used  to  measure  voltage.  Ohe  voltmeter  with  a  range  of  0  to 
50.000  volts,  was  incorporated  as  part  of  the  power  supply  and  was  read  during  all  heat 
transfer  runs.  The  second  voltmeter  an  RCA  Voltohmyst,. utilized  a  high-voltage  probe, 
and  drew,  undesirably  high  currents  during  operation;  these  currents  could  be  accounted 
for.  however,  and  calibrations  were  set  up  to  correct  for  them.  The  third  voltmeter,  a  . 
Sensitive  Research  Instmment  Corporation  electrostatic  voltmeter  with  a  range  of 
.0-40,000  volts,  was  used,  during  precise  voltage,  current,  and  pressure  measurements 
of  various  corona  discharges.  Independent  of  wWch  voltmeter  was  used  in  the  tests,  it 
was  foiuid  that  the  current  was  the  most  sensitive  parameter  to  control  and,  as  a  conse¬ 
quence,  all  settings  were  based  upon  current.  It  was  often  noted  that  large  variations  in 
current  occurred  even  though  the  changes  in  voltage  were  barely  discemable. 

Fot  the  tests  with  dielectrophoresis,  an  A.C.  power  supply  was  used.  The  supply  con¬ 
sisted  of  a  neon  transformer  connected  to  a  Variac  to  provide  a  voltage  range  of  zero  to 
20,000  volts.  Several  milllamperes  of  current  could  be  drawn  from  this  source. 
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Electrodes 


Several  types  of  electrodes  that  were  fabricated  and  tested  are  shown  in  Figure  7,  The 
electrode  at  the  left  consists  of  a  series  of  No.  7  needles  soldered  precisely  at  1/2-inch 
intervals  along  a  1/8-lnch  brass  rod  11  Inches  long.  The  tips  of  the  needles  were  aligned 
carefully  in  a  jig  so  that  distortion  in  the  interferometer  due  to  inaccuracies  of  needle 
arrangement  would  be  at  a  minimum.  The  needles  in  this  case  served  as  tte  corona  dis¬ 
charge  points.  The  nondischarging  electrode  at  the  front  Is  made  of  1/8-inch  brass  rod 
shaped  into  a  12-inch  oval.  The  electrode  In  the  rear  of  the  picture  is  a  grid  of  fine  wires 
0.010  Inch  in  diameter,  12-inches  long,  and  spaced  at  one-inch  intej^als.  The  wires  act 
as  the  corona  discharge  source.  The  electrode  at  the  right  is  an  unheated,  flat,  6  by  10- 
inch  aluminum  plate  that  was  used  as  a  nondischarging  electjrode  to  provide  an  approxi¬ 
mately  parallel  field. 

The  installation  of  an  electrode  typical  of  those  used  In  these  tests  js  shown  iii  Figure 
5,  where  the  heated  6-by-lO-Inch  aluminum  plate  is  mounted  in  the  interferometer.  A  • 
single-wire  electrode  can  be  seen  mounted  at  a  distance  of  2  centimeters  in  front  of  the 
plate.  This  single-wire  electrode  is  typical  of  the  type  used  primarily  in  these  tests.  It 
is  a  single  stainless-steel  wire  kept,  under  tension  to  keep  it  'straight  end  electrically  In-  • 
sulated  from  the  mounting  attachments.  A  small  stranded-wire  lead  extends  from  the  ' 
corona  wire  out  to  the  point  of  attachment  of  the  high-voltage  lead.  The  electrode  is 
mounted  in  a  traversing  mechanism  that  allows  movement  in  three  perpendicular  co¬ 
ordinate  directions.  Thus  the  corona  wire  may  be.  moved  up  and  down  or  Its  distance 
from  the  plate  varied  as  required.  The  wire  shown  In  the  fi^te  ts  12  inches  long  and 
0.004  inch  in  diameter.  ...  .  '  . 

Both  the  plate  and  the  electrode  were  electrically  Insulated  from  th^  remainder  of  the  ' 
equipment.  The  plate  was  maintained  at  ground  potential  by  the  ground  Vetum  lead  to  the 
power  supply.  All  the  sharp  edges  of  the  plate  and  the  electrodes  were  rounded  or  covered 
with  cement  or  corona  dope  to  minimize  undesirable  localized  corona  discharge,  and  the 
surface  of  the  plate  was  burnished  and  polished,  clean  and  free  of  sharp  scratches.  In 
spite  of  the  precautions,  however,some  localized  discharge  occurred  both  from  spots  on 
the  plate  and  from  distinct  points  on  the  various  electrodes- at  very  high  values 'of  applied 
voltage.  Fortunately,  these  localized  discharges  appeared  to  occur  only  at  voltages  and 
Currents  beyond  the  range  of  Iitterest  in  the  study. 

Pressure  Survey  Equipment 

To  interpret  the  phenomena  observed  with  the  interferometer,  it  became  iTecessary  to 
study  the  pressure  distribution  between  various  electrodes  and  the  [flat  plate.  Review  of 
data  obtained  by  Chatlock,  Harney,  Stuetzer,  and  Robinson  (Refs,  5,  6.  7.-  and  8.  respec¬ 
tively)  revealed  that  a  sharp  pressure  peak  could  be  expected  directly  beneath  the  fine 
wire  during  a  corona  discharge.  None  of  these  investigators,  hpwever,  had  published  .any 
data  for  the  case  of  fine  wires  placed  parallel  to  a  flat  plate.  , 

To  determine  the  pressure  distribution;  tests  were  run  using'an  arrangement  similar 
to  that  shown  in  Figure  5,  but  the  heated  plate  was  replaced  by  a  pressure  plate  equipped 
with  a  series  of  pressure  taps.  Pressures  were  recorded  during  the  corona  discharge, 
but  no  interferometer  or  heat  transfer  data  were  taken  during  the  pressure  surveys; 

Two  different  test  configurations  , to  measure  pressure  variation  were  set  up.  In  the 
first  configuration,  a  manometer  board  was  tilted  to  a  10°  angle.  Flood  lights  were'  used 
for  sharp  definition  of  the  meniscus  for  all  the  tubes,  as  shown  in  Figure  8.  To  measure 
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the  pressures,  a  series  of  holes  were  drilled  in  a  7  by  12  inch  aluminum  plate  at  0.2-inch 
intervals  with  a  No.  54  drill  along  the  horizontal  and  vertical  centerlines  of  the  plate.  One 
hole  was  drilled  at  the  center  of  the  plate  and  the  remainder  were  spaced  outwardly  to  a 
distance  of  two  inches  from  the  center.  The  back  of  the  plate  was  counterdrilled  so  that 
1/16  inch  metal  tubing  could  be  inserted.  All  tubes  were  sealed  into  the  plate  with  a  Glyptol 
varnish.  The  pressure  measuring  plate  is  shown  at  the  top  of  the  manometer  board  in  Fig¬ 
ure  8. 

Plastic  tubing  was  used  to  connect  the  pressure  taps  to  the  manometer  board.  All  leads 
were  made  as  short  as  possible,  and  all  connections  were  made  as  tight  as  possible  to 
minimize  the  possibility  of  leakage.  The  entire  surface  of  the  plate  was  polished.  All  nicks, 
scratches,  and  rough  edges  were  smoothed  out,  and  corona  dope  was  applied  at  all  positions 
of  potential  local  discharge  away  from  the  center  of  the  plate. 

To  obtain  useful  data  on  the  very  small  pressure  differences,  special  procedures  had  to 
be  employed  in  using  the  mancMneter  board.  Room  temperature  changes  had  to  be  watched 
carefully  to  minimize  drift  effects.  A  Polaroid  camera  was  mounted  above  and  perpendic-  • 
ular  to  the  manometer  board  to  minimize  distortion.  Polaroid  pictures  were  taken  of  each 
test  sequence  and  read  by  means  of  a  Gaertner  Comparator.  Reference  pictures  were  taken 
before  and  after  each  test  sequence.  By  comparing  the  position  of  the  meniscus  on  each 
data  run  with  the  corresponding  reference  values,  it  was  possible  to  read  values  of  pres¬ 
sure  to  less  than  0.001  inches  of  water.  Repeatability  was  good.  The  simultaneous  readings 
of  pressure  across  the  face  of  the  plate  gave  a  good  indication  of  the  form  of  the  pressure 
distribution. 

The  second  technique  of  measuring,  pressure  utilized  the  same  electrodes  and  pressure 
plate,  but  the  pressures  were  read  hy  means  of  a  Chattock  micromanometer,.  T.he  photo¬ 
graph  of  this  pressure  test  set-up.  Figure  9,  shows  the  high-voltage  power  supply,  elec¬ 
trostatic  voltmeter,  and  ammeter,  as  well  as  the  pressure  plate  and  manometer;  Since 
the  Chattock  micromanometer,  shown  in  Figure  10,  provides  only  one  pressure  reading 
at  a  time,  it  was  necessary  to  change  the  connection  between  the  pressure  taps  and  the 
manometer  in  sequence.  The  plastic  tubing  lead  was  shortened  as  much  as  possible  to 
minimize  the  lag  in  the  system  and  the  lines  made  as  clear  as  possible.  Readings  to  de-  • 
termine  the  lag  time  indicated  that  a  stable  value  was  achlevjsd  well  within  one  minute. 

The  Chattock  manometer  uses  a  two-fluid  “U”  tube.  The  fluids  used  in  this  base  were 
kerosene  and  distilled  water.  A  Gaertner  Comparator,  which  was  used  to  read  the  .level 
of  the  meniscus,  increased  the  sensitivity  of  the  instrument  and  facilitated  reading.  The 
Chattock  manometer,  however,  requires  considerable  adjustment  and  care  in  Operation  if 
consistent  data  are  to  be  obtained.  Rown  temperature  variations,  outside  winds,  and  gusts 
readily  affected  not  only  the  individual  readings  but  the  zero  shifts,  and  variations  in  hu¬ 
midity  affected  the  corona  discharge.  Consequently  pressure  tests  were  run -only  on  quiet 
days  with  the  room  heat  turned  off,  and  on  days  when  the  relative  humidity  was  not  high. 
Figure  11  illustrates  a  typical  set-up  with  the  single  wire  in  place  a  short  distance  in 
front  of  the  pressure  plate.  .  • 

Experimental  Procedures 

Because  electrostatic  phenomena  produce  forces  that  are  small  in  magnitude,  it  was 
necessary  to  watch  for  all  extraneous  disturbances  that  could  occur  in  the  test  environ¬ 
ment  as  well  as  the  accuracy  of  the  instrumentation.  The  difficulties  with  the  temperature 
variations,  gusty  weather,  and  hurnidity  have  been  indicated  previously.  To  alleviate  the 
presence  of  small  random  air  motions  in  the  vicinity  of  the  test  section  of  the  interfero¬ 
meter,  an  enclosure  was  placed  about  this  area.  This  enclosure  was  located  so  as  to  not 
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materially  affect  the  free -convection  flow  of  air  near  the  plate,  the  flow  of  air  from  the 
plate,  and  the  flow  of  air  from  the  corona  discharge.  Use  of  the  enclosure  resulted  in 
more  stable  test  patterns. 

An  interferometer  provides  data  on  the  isothermal  lines  in  the  fluid  around  a  heated 
body  by  sensing  the  average  density  changes  in  the  working  section.  A  free-convection 
test  of  a  flat  plate  in  still  air  will  reveal  a  series  of  fringes  in  front  of  the  plate.  By 
measuring  the  spacing  of  the  fringes  and  calculating  the  temperatures  existing  in  each 
fringe,  it  is  possible  to  determine  the  temperature  profile  through  the  thermal  boundary 
layer  as  well  as  the  local  heat  transfer  coefficients.  A  typical  interferometer  picture. 
Figure  12,  shows  an  end  view  of  the  test  specimen  with  the  flat  heated  surface  at  the  right. 
The  fringes  are  evident  in  front  of  the  plate.  In  this  case,  the  plate  temperature  was  92*F 
and  the  room  temperature  was  64.1*F.  The  “T”-shaped  rod  at  the  bottom  provides  a  ref-^ 
erence  length  for  data  reduction.  The  end  view  of  the  electrode  support  rod  can  be  seen  at 
the  left.  side.  The  fine  corona  wire  at  the  end  of  the  rod  cannot  be  discerned  in  the  picture. 
No  voltage  was  applied  at.  the  time  this  interferogram  was  taken.  When  a  voltage  is  applied 
between  the  fine  wire  and  the  plate  and  changes  to  the  thermal  boundary  layer  occur,  the 
fringe  patterns  are  distorted.  In  general,  if  the  fringes  are  brought  closer  together  so  that 
the  thermal  boundary  layer  is  thinner,  heat  transfer  increases.  If  the  boundary  layer  is 
pulled  outwardly,  the  heat  transfer  rate  decreases.  Thus,  when  viewing  the  data  presented 
in  subsequent  [^otographs,  these  basic  points  should  be  borne  in  mind  as  an  aid  in  inter¬ 
preting  the  phenomena. 

The  accuracy  of  the  data  depends  on  both  the  environmental  conditions  and  the  instru¬ 
mentation  used.  Temperatures  measured  with  the  thermocouples  and  the  thermometers 
could  be  read  to  within  0.1“F.  The  interferometer  fringe  patterns  contain  errors  due  to 
the  variation  in  vlenaity  In  the  horizontal  (Mreetlon  aking  the  {gate.  Ihis  eirror  can  be 
compensated  for  and  a  compensation  is  included  in  the  data  reduction  method.  A  second 
error  is  introduced  by  the  end  effects  of  the  finite  plate  width.  The  largest  inaccuracy 
conU'B  fa  actual  n^sauainninJi  of  the  apactfigi  Itw’  friagiiA.  Alttxiughthe  Caermer 
Comparator  can  measure  the  spacing  to  O.OOOl  inch,  determining  the  centerline  of  .the 
fringes  introduces  far  greater  uncertainties.  When  the  fringes  are  diffused  and  spread 
far  apart,  a  sizeable  unceitalnty  can  result.  Conversely,  an  uncertainty  is  also  introduced 
when  the  fringes  are  compacted,  one  upon  the  other,  because  the  Individual  lines  are  dif¬ 
ficult  to  locate. 

~  The.  voltage  and  current  data  are  also  subject  to  Inaccuracies.  As  mentioned  previously, 
the  voltage  Used  In  corona  discharge  phenomena  Is  a  relatively,  insensitive  parameter  in 
the  region  of  the  initial  flow  of  corona  current.  In  order  to  initiate  the  corona  discharge, 
the  voltage  control  knob  of  the  power  supply  was  turned  slowly  until  a  change  in  current  . 
could  be  observed,  but  the  voltage  reading  was  of  little  assistance  in  detemining  the 
starting  point.  Many  calibrations  were  made  with  the  three  voltmeters  described  previously. 
The  precision  electrostatic  voltmeter,  which  was  used  as  a  standard  to  calibrate  the  others, 
was  available  for  use  during  Only  a  portion  of  the  work;  as  a  consequence,  recourse  was 
made  to  the  less  accurate  instruments. 

Current  measurements  made  with  the  Greibach  Microammeter  are  considered  to  be  • 
accurate.  Unfortunately,  corona  discharge  currents  are  very  erratic.  The  fundamental 
nature  of  the  discharge  currents  with  a  positive  point  consists  of  a  series  of  random 
current  pulses  accompanying  electron  avalanches  towanl  the  positive  point.  Slight 
variations  in  atmospheric  conditions  in  the  yicinity  of  the  wire  could  also  cause  variations 
.  in  the  corona  current.  Consequently,  the  pointer  of  the  microammeter  moved  continually 
and  erratically  during  all  readings.  These  same  variations  have  been  noted  by  other 
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investigators  (Refs.  6,  8).  Several  readings  were  usually  taken,  but  drifting  would  occur 
to  either  the  high  or  low  side,  and  the  averages  were  not  truly  representative.  Efforts 
were  made  In  all  tests  to  take  readings  at  the  most  stable  current  value  in  the  range  under 
consideration.  Because  of  this  erratic  behavior,  it  was  difficult  to  assign  a  value  to  the 
accuracy  of  the  current  data  even  though  it  is  a  most  sensitive  .parameter  in  this  investi¬ 
gation.  Prom  study  of  many  data  points,  it  is  believed  that  the  current  values  are  within 
5  per  cent  in  the  small  current  ranges  (less  than  10  4  a)  and  considerably  better  than  that 
in  the  higher  ranges. 

The  pressure  data  ctmtain  several  possible  sources  of  inaccuracies.  With  both  mano¬ 
meters,  the  possibility  existed  of  the  meniscus  sticking  to  the  wall  of  the  tubes  and  thus 
introducing  errors.  If  light  were  directed  continuously  on  either  manometer,  all  readings 
became  meaningless;  volumetric  changes  and  distortions  of  the  apparatus  occurred  from  * 
exposure  to  direct  light  for  more  than  an  instant.  Consequently,  when  working  with  the 
Inclined  manometer,  the  lights  were  flashed  on  momentarily  while  taking  a  picture  of 
the  tubes. 

The  Chattock  manometer  was  particularly  sensitive  to  random  external  disturbances. . 

As  Indicated  before,  changes  in  temperature,  drafts,  outside  winds,  and  aircraft  in  flight 
all  could  cause  erratic  readings.  In  addition  to  these,  even  such  disturbances  as  doors 
being  closed,  cabinets  opened,  and  other  tests  being  conducted  in  remote  areas  could 
cause  severe  transients  in  the  pressure  readings.  An  ideal  location  In  which  to  perform 
these  tests  would  have  been  a  room  with  air  conditioning,  controlled  humidity,  extremely 
rugged  walls,  and  heavy  thermal  Insulation. 

By  setting  the  pressure  equipment  up  in  the' same- room  with  tte  interferometer,  dis- 
luibances  due  to  fhe  movement  of  peujile  was.  buhstanllally  eliminated.  Testing  was  accom¬ 
plished  on  what  appeared  to  be  very  quiet  and  calm  days.  Even  under  the  best  conditions, 
a  troublesome  drift  of  the  referees  zero  occurred.  Many,  runs  were  obtained  under  good 
test  conditions  In  which  the  zero  shift  was  within  O.QOOl  inch  of  water,  and  these  data 
show  a  consistent  pattern  with  little  shifting  and -distortion.  Unfortunately,  .zero  shifts 
occurred  during  some  of  the  runs  that  consumed  several  hours,  during  which  time  siftal! 
atmospheric  changes  may  have  occurred  which  were  not  readily  detectable.  Data  obtained 
under  these  conditions  were  erratic  and  did  not  correlate;  weU.  In  spite  of  the'  difficulties, 
we  believe  that  the  data  as  presented  .in  the  next  section  are  reliable,  ,  .  • 

By  far  the  most  troublesome  problem  in  the  whole  experiment  was  the  difficulty  in 
repeating  tests  and  obtaining  data  with  a  high  degree  of  precision  (within  oiie  percent). 

Some  of  the  factors  causing  this  have  been  mentioned  before,  and  these  are  of  consider- 
able  Importance.  Fundamental  to  this  prtAlem,  however,  is  the  erratic  character  of  the 
corona  discharge  Itself.  Variations  of  the  corona  due  to  small  changes  in  humidity,  local  ' 
ionization,  and  random  convective  currents  are  very  difficult  to  control.  The  average 
mobility  of  the  ions  is  affected  considerably  by  the  amount  of  moisture  in  the  air.  Since 
pressure  rise  generally  varies  inversely  with  Ion  mobility,  pressure  rise  due  to  corona 
discharge  can  be  expected  to  vary  with  humidity.  Since  clear-cut  corrections  for  the 
influence  of  humidity  on  pressure  rl^  were  not  available,  the  pressure  rise  data  were 
not  corrected  for  humidity  effects.  The  basic  interactions  between  corona  discharge  and 
ho-s(  Irmutcr,  (mvever,  ocourzvJ  at  all  of  humiditf,  from  very 

dry  cold  days  to  foggy  days  with  the  humidity  close  to  100  percent. 
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RESULTS  OF  THE  EXPERIMENTAL  INVESTIGATION 
Method  of  Data  Interpretation 

The  heat  transfer  coefficients  were  obtained  from  the  interferograms  by  measuring 
the  spacing  of  the  fringes  and  using  the  calculation  procedures  presented  In  Appendix  1. 
Readings  were  taken  from  at  least  ten  separate  locations  along  each  plate.  The  data  are 
presented  in  plots  to  show  heat  transfer  coefficients  along  the  surface  of  the  plate,  and 
the  effects  of  varying  current  and  wire  spacing.  Voltage -current  graphs  are  presented 
to  illustrate  the  usual  characteristics  of  corona  discharge. 

Interferometer  Data 

The  results  of  the  combined  electric  field-convection  tests  are  shown  in  the  accompany¬ 
ing  photographs  and  correspcxiding  graphs.  Only  the  most  significant  results  are  presented 
and  interpreted.  The  first  sequence  of  tests  Incorporates  the  heated  plate  in  a  vertical  po¬ 
sition,  with  a  single  stainless-steel  electrode  0.004  inch  iii  diameter  located  2  cm  away 
from  the  plate.  This  configuration  provides  an  intense  electrical  field  near  the  wlrej  and 
a  highly  nonuniform  electrical  field  in  the  region  in  front  of  the  plate.  The  wire  is  hori¬ 
zontal  and  at  the  centerline  of  the  plate.  The  potential  applied  between  the  wire  and  the 
plate  is  slowly  raised  from  zero  up  to  values  sufficient  to^give  currents  of  abproximatelv 
150pa.  .  ■  . 

Figure  13  presents  a  sequence  of  pictures  as  seen  with  the  Interferometer  as  the  volt¬ 
age  is  increased.  View  A  is  the  reference  photograph  with  no  applied  field;  plate  tem¬ 
perature  is  119'’F  and  room  temperature  is  69®F.  In  this  exploratory  run.  .no  attempt  was 
made  to  hold  the  plate  temperature  constant  at  high  corona  current  values.  At  high  values 
of  corona  current,  plate  temperature  dropped  to  approximately  lOO^F  with  a  constant  heater 
current.  The  fringe  pattern  shown  in  this  picture  remains  unchanged  as  the  potential  Is 
applied  as  long  as  the  current  remains  essentially  Zero.  As  soon  as  a  slight  current  begins 
to  flow.  Some  changes  in  the  pattern  occur.  This  behavior  is  characteristic  of  all  tests 
conducted  with  the  corona  discharge.  The  rods  holding  the  fine  wires  are  visible  at  the 
left-  The  “T”  located  at  the  bottom  of  the  picture  provides  a  reference  length. 

The  initial  disturbance  of  the  thermal  boundary  layer,  shown  in  View  B,  occurs  just  at 
the  onset  of  corona  current.  Subsequent  pictures,  in  sequence,  show  the  marked  distortion 
of  the  Isothermal  lines  as  the  current  is  increased.  The  patterns  shown  are  fairly  steady 
for  a  giyen  stabilized  current  except  for  Views  D  and  E.  which  are  slightly  unsteady.  The 
isothermal  lines  are  extremely  steady  for  currents  in  excess  of  10  u  a.  The  pictures  reveal 
that,  depending  upon  the  value  of  current,  regions’ of  Increased  or  decreased  heat  transfer 
exist.  The  data  on  heat  transfer  coefficients,  hx,  in  this  sequence  are  considered  quite  ac¬ 
curate  for  values  of  current  less  than  10  jua.  Above  this  value,  sufficient  plate  cooling  oc¬ 
curred  to  introduce  some  error  in  hx.  Consequently,  no  hx  data  are.  presented  fox  currents 
greater  than  10  p  a.  The  values  of  the  heat  transfer  coefficients  corresponding  to  Views  A 
through  E  are  plotted  in  Figure  14.  The  heat  transfer  coefficient  at  various  locations  along 
the  plate  is  plotted  for  several  values  of  corona  discharge  Current.  The  abscissa,  x,  rep¬ 
resents  the  distance  in  the  vertical  direction  along  the  plate  surface;  x  =  0  represents  the 
centerline  of  the  plate,  three  inches  from  the  bottom  edge;  distances  to  the  right  of  x  =  0 
represent  the  lower  part  of  the  plate,  and  distances  to  the  left  represent  the  upper  part  of 
the  plate.  A  marked  change  in  the  magnitude  and  distribution  of  hx  results  from  the  appli¬ 
ed  electrical  field.  For  small  currents,  the  variation  along  the  plate  is  extremely  large, 
while  for  higher  currents,  the  values  of  hx  tend  to  become  more  symmetric  about- the 
centerline. 
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Figure  15  presents  the  results  of  tests  using  the  same  configuration  as  in  Figure  13  but 
the  temperature  difference  between  plate  and  room  air  is  maintained  constant  throughout 
the  test  sequence.  Stabilized  test  points  were  obtained  and  multiple  records  were  made  of 
the  phenomena.  Accurate  data  were  obtained  up  to  high  values  of  current.  View  A  shows 
the  reference  pattern  with  no  applied  field.  Changes  in  the  fringes  are  evident  as  the  cur¬ 
rent  is  increased.  At  the  higher  values  of  current,  the  thermal  boundary  layer  becomes 
extremely  thin;  it  can  be  expected,  therefore,  that  the  heat  transfer  rate  has  increased 
greatly.  Unfortunately,  this  extreme  thinness  of  the  boundary  layer  makes  It  difficult,  if 
not  impossible,  to  determine  the  local  values  of  heat  transfer  coefficient.  The  heat  trans¬ 
fer  coefficients  corresponding  to  this  sequence  of  photographs  are  plotted  in  Figure  16. 

The  reference  curve  with  no  field  applied  can  be  obtained  with  considerable  accuracy. 
Values  of  heat  transfer  coefficient  for  the  higher  values  of  current  indicate  a  larger  de¬ 
gree  of  scatter,  and  could  not  be  determined  for  the  center  of  the  plate. 

The  next  series  of  tests,  shown  in  Figure  17,  was  conducted  with  the  electrode  moved 
away  from  the  plate  to  a  distance  of  5  cm.  The  temperature  difference  was  maintained  at 
34°F.  The  thermal  boundary  layer  is  squeezed.down  as  the  field  is  increased.  The  main 
difference  between  this  case  and  the  previous  one  at  corresponding  high  currents  is  that 
the  thermal  layer  is  thinned  more  uniformly  across  the  plate.  Corresponding  values  of 
heat  transfer  coefficient  are  plotted  in  Figure  18.  Difficulties  were  again  encountered  in 
reducing  the  data  at  the  higher  current  values. 

A  test  sequence  with  the  electrode  moved  out  to  io  cm’ is  shown  in  Figure  19.  The  high¬ 
er  currents  were  not  reached  in  this  case  because  the  required  voltage  became  excessively 
high.  At  high  voltages,  local  points  of  corona  discharge  occurred  on  the  plate  and  on  dis¬ 
crete  points  of  the  wire.  Such  intense  local  discharges  gave  erroneous  impressions  of  the 
value  of  the  current  driving  the  corona  wind.  Heat  transfer  coefficients  are’ shown  In  Fig¬ 
ure  20.  ' 

In  the  next  sequence  of  tests,  the  heated  plate  was  placed  with  the  flat  heated  surface 
faced  downward  and  the  0.004  inch  wire  electrode  2  cm  away  from  the  plate  at  the  center- 
line.  The  temperature  difference  was  held  constant  at  27°F.  The  influences  ■on  the  isotherms 
as  the  Current  begins  to  flow  for  the  range  from  0  to  4.55  fia  only  are  Indicated  In  Figure 
21.  The  gradual  squeezing  of  the  bpundary  layer  is  clearly  evident.  Other  carefully  con¬ 
trolled  experiments  indicated  that  the  patterns  shown  is  very  symmetrical  about  the  center- 
line  for  all  values  of  current.  The  corresponding  heat  transfer  coefficients  are  plotted  in 
Figure  22.  The  method  of  presenting  data  four  the  horizontal  plate  Is  similar  to  that  used 
for  the  vertical  plate.  Since  the  configiira  tion  was  symmetrical  about  the  centerline  of  the 
plate,  however,  data  for  only  one  half  of  the  plate  are  plotted,  from  the  centerline  out  to  • 
the  edge.  The  gradual  ctange  in  magnitude  and  distribtuion  Is  clearly  evident.  The  varia¬ 
tion  of  hx  becomes  bell-shaped  as  the  current  is  applied.  At  the  edge  of  the  plate,  where 
x  =  3,  the  end  effects  cause  a  localized  increase.  This  characteristic- increase  is  eyident 
in  the  curve  for  free  convection  with  no  field  applied.  Ir  is  interesting  to  notice  that  the 
local  values  of  hx  can  be  changed  several  fold  with  currents  of  a  few  microamperes. 

Figure  .23  presents  the  results  of  another  test  sequence  in  which  the  same  configuration 
was  used,  but  the  temperature  difference  was  maintained  at  27“F  and  the  current  range 
covers  0  pa  to  159  ua.  The  thinning  of  the  thermal  boundary  layer  Is  quite  evident.  At  the 
high  values  of  current,  the  patterns  resemble  those,  found  with  the  vertical  plate.  The 
heat  transfer  coefficients  are  plotted  in  Figure  24.  The  values  at  high  currents  and  near 
the  center  tend  to  scatter  as  before. 
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In  the  next  sequence  of  tests,  the  same  configuration  of  horizontal  plate  is  used,  but  the 
wire  is  moved  to  a  distance  of  5  cm  and  the  temperature  difference  is  maintained  at  34°F. 
The  results  of  the  tests,  shown  in  Figure  25,  again  indicate  a  thinning  down  of  the  bound¬ 
ary  layer,  but  the  shape  is  somewhat  flatter  than  that  with  the  electrode  at  2  cm.  The 
heat  transfer  coefficients  plotted  in  Figure  26  show  the  changes  in  magnitude  and  distri¬ 
bution.  Once  again,  at  higher  current  values  the  data  becomes  scattered  due  to  the  diffi¬ 
culty  of  reading  the  interferograms.  Large  changes  in  the  heat  transfer  coefficient  again 
result  from  a  small  current  flow. 

Pressure  Survey  Data 


Tests  of  the  pressure  rise  obtained  with  a  corona  wind  were  conducted  to  provide  data 
on  the  magnitudes  of  the  pressures  as  well  as  the  variation  of  the  pressure  along  the  plate. 
The  test  configuration  consisted  of  locating  a  single  wire  electrode  parallel  to  a  7  by  12- 
inch  flat  plate  at  the  centerline.  Wires  measuring  0.0005,  0.004,  and  0,010  inch  in  diameter 
were  tested  at  various  distances  from  the  plate.  The  wire  was  oriented  so  that  it  ran,  par¬ 
allel  to  the  7-inch  dimension  of  the  pressure  plate.  The  pressure  rise  obtained  at  the  cen¬ 
terline  of  the  plate  with  the  different  size  electrodes  is  shown  in  Figures  27,  28,  and  29. 

All  plots  tndicate  an  approximately  linear  preBSure-current  characteristic,  which  is  in 
accordance  with  the  theories  of  corona  wind.  Some  zero  shift  occurred.  A  typical  zero 
shift  can  be  seen  in  the  2-cm  case  of  Figure  29.  A  small  influence  of  spacing  of  the  wire 
can  be  noted,  but  trends  due  to  wire  spacing  are  not  consistent  as  the  wire  size  is  changed. 
The  magnitudes  of  the  pressure  for  a  given  current,  likewise,  do  not  seem  to  be  markedly 
influenced  by  wire  size.  No  clear-cut  trends  of  the  effects  of  wire  size  or  wire  spacing 
can  be  determined  from  the  data.  AH  data  presented  in  these  three  figures  were  taken 
with  the  Chattock  micromanometer. 

The  next  sequence  of  tests  were  conducted  to  investigate  pressure  distribution.  Data 
were  taken  with  the  Chattock  micromanometer,  A  0.004-inch  wire  electrode  was  used  at 
three  wire  spacings,  2  cm,  4  cm,  and  6  cm.  Results  of  these  tests  are  plotted  in  Figures 
30,  31  and  32.  The  magnitudes  of  the  centerline  pressures  at  a  given  current  do  not  vary 
greatly,  as  had  teen  indicated  in  the  previous  test  data.  The  shape  of  the  curves  at  all 
spacings  is  similar.  Since  testing  had  shown  that  the  pressure  distributions  were  symmet¬ 
rical,  only  one-half  of  the  distribution  is  shown.  Because  the  corona  discharge  does  not 
necessarily  occur  in  a  line  perpendicular  from  the  wire  to  the  plane,  some  shift  in  the 
peak  pressure  position  can  be  expected.  Study  of  the  curves,  however,  reveals  that  the 
centerline  of  the  corona  wind  actually  did  move  away  from  the  geometric  centerline  of 
the  plate  in  some  cases. 

Figure  33  presents  the  results  of  another  series  of  tests  at  a  2-em  spacing  to  show  the 
characteristic  variability  of  the  corona  phenomenon.  Data  taken  at  the  lower  three  cur¬ 
rents  were  strongly  affected  by  random  zero  shifts.  Although  the  curves  seem  to  be  of 
reasonable  proportion,  a  subsequent  curve  shows  that  these  data  are  unusable., 

,  Typical  pressure  distribtuions  obtained  with  the  inclined  manometer  are  presented  in 
Figure  34  to  show  the  general  shape  of  the  pressure  distribution  curve.  The  absolute 
values  of  these  data  are  not  as  accurate  as  those  obtained  with  the  Chattock  manometer, 
but  the  curve  has  the  same  tell  shape. 

Current-Voltage  Relations 

The  variation  of  current  with  a*^ied  voltage  was  meaRiired  often  during  XIk.  niUirnii.i  uT  ' 
the  testing.  Since  it  proved  difficult  to  relate  the  voltage  to  either  test  technique  or  data 
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interpretation,  only  one  representative  voltage  curve  will  be  presented.  Figure  35  pre¬ 
sents  the  results  for  a  corona  discharge  with  a  0.004-inch  wire  located  at  the  centerline 
of  a  7  X  12-inch  flat  plate  and  running  in  the  direction  of  the  7-lnch  dimension.  Data  are 
presented  for  various  wire  spacinge.  It  is  believed  that  the  results  are  typical  of  those  . 
found  with  corona  discharge. 

Discussion  of -Test  Results 

The  purpose  of  this  section  is  to  desciribe  in  some  detail  the  phenomena  observed  in  the 
foregoing  sequence  of  tests  on  the  Interferometer. 

An  unusual  movement  of  the  thermal  boundary  layer  is  shown  in  Figure  13.  This  phe¬ 
nomenon  occurred  in  all  runs  with  the  vertical  plate,  and  with  several  electrode  configu¬ 
rations.  At  the  initial  flow  of  current,  the  isotherms  are  “pulled-out"  from  the  heated 
plate.  At  soriiewhat  greater  currents  this  pull-out  becomes  larger.  As  the  current  was 
raised  to  even  higher  values,  the  pull-out  disappeared.  The  cause  of  the  pull-out  was  in¬ 
vestigated  in  considerable  detail.  Pressure  distribution  measurements  at  low  currents 
failed  to  provide  any  accurate  data  on  the  pressure  field  in  the  region  of  pull-out.  The  pull¬ 
out  shape  appeared  to  move  opposite  to  the  direction  expected  front  corona  wind  action. 

Any  dipole  forces  acting  in  the  air  should  also  have  acted  in  a  direction  opposite  to  the 
puU-out.  A  possible  explanation  of  the  piill-out  phenomenon  lay  in  the  action  of  two  im¬ 
pinging  fluid  streams.  With  a  heated  vertical  plate,  a  convective  stream  moves  upward 
along  the  plaite.  If  the  corona  wind  strikes  the  plate,  it  can  be  expected  that  the  air  stream 
will  split  and  try  to  move  upward  and  downwarf.  The  corona  stream  moving  downward  will 
buck  the  upward  streaming  of  the  thermal  boundary  layer.  It  appears  possible,  therefore, 
that  as  the  two  fluid  streams  meet,  the  temperature  field  would  change  in  an  unanticipated 
manner.  The  possibility  that  the  pull-out  phenomenon  was  due  to  an  interaction  of  the  cor¬ 
ona  wind  and  the  thermally  heated  stream  was  explored  in  the  sequence  of  tests  shown  in 
Figures  21,  23,  and  25.  The  corona  wire  was  placed  exactly  at  the  center  of  an  accurately 
ievel^  horirontai  plate  at  a  short  distance  below  the  plate.  Upon  applicatton  xA  the  field,  . 
the  center  of  the  corona  wind  stream  would  be  located  at  the  point  of  zero  longitudinal 
velocity  of  the  thermal  boundary  layer.  Thus,  any  corona  stream  striking  the  plate  would 
split  into  two  streams,  each  of  which  would  move  with  the  streaming  in  the  thermal  bound¬ 
ary  layer.  If  no  pull-out  occurred,  then  the  pull-out  phenomenon  must  have  been  associated 
with  the  fluid  flow  impingement.  Conversely,  if  pull-out  did  occur,  then  the  phenomenon 
requKed  fiome  «hd  this,  bk  %  mem  ^ps  dl  4m&¥h(%lor>. 

The  tests  with  the  horizontal  plate  were  setup  with  the  hot  side  up  and  the  hot  side  down. 
The  convective  cells  developed  in  free  convection  with  the  hot  face  up  made  it  impossible 
to  interpret  the  results.  Tests  with  the  hot  side  down  and  the  corona  wire  at  several  spac- 
ings,  however,  were  most  successful,  as  shown  in  Figures  21,  23,  and  25.  As  field  is 
applied  and  the  current  flows,  no  pull-out  is  evident  with  either  the  2  cm  or  5  cm  wire 
spacing.  The  experiments  were  repeated  several  times  to  verify  results,  and  extreme 
caution  was  used  so  that  extraneous  disturbances  did  not  introduce  error.  Since  the  pull¬ 
out  did  not  occur,  it  was  felt  that  the  pull-out  phenomenon  observed  with  the  venical 
plate  tests  was  due  to  the  impingement  of  fluid  streams,  and  not  to  an  unexplained  elec¬ 
trostatic  effect. 
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ANALYTICAL  STUDY  OF  THE  EFFECT  OF  CORONA  DISCHARGE 
ON  FREE-CONVECTION  HEAT  TRANSFER 


The  influence  of  the  applied  electric  field  on  the  heat  transfer  of  the  heated  flat  plate 
was  analyzed  in  an  attempt  to  learn  possible  trends  that  might  occur  and  to  aid  in  obtain¬ 
ing  a  better  understanding  of  the  phenomena.  The  analytical  work  is  broken  down  into  four 
phases:  electrical  field  effects,  heat  transfer  theory,  the  combined  problem^and  an  approx¬ 
imate  solution  to  the  horizontal  plate  case. 

ELECTRIC  field-corona  DISCHARGE  RELATIONS 

In  analyzing  the  corona  discharge  phenomena,  it  is  essential  to  obtain  relationships 
covering  the  field  strength,  current  distribution,  current-voltage  relationships,  and  the 
induced  pressure  rise.  The  theory  is  based  upon  the  recognition  that  a  definite  charge. 

Pc  exists.  Two  simple  illustrations  of  solutions  will  be  presented:  a  one-dimensional 
case,  and  a  cylindrical-electrode  configuration. 

The  electric, field  equations  (from  Ref.  1)  are: 

V  •  E  =  — ^ 

< 

T  =  oc  .  E  +  7^ 

F  =  E 

Pressure  Rise  with  Space  Charge 

_  Ih'a  gas  where  the  total  current  carriers  are  associated  with  the  ions,  J  *  pc  v^,  where 
Vj  is  the  total  velocity  of  the  moving  charge.  It  is  assumed  that  ions  of. only  one  sign  exist 
away  from  the  immediate  vicinity  of  the  corona  point,  and  the  mobility  of  these  ions  is 
approximately  constant.  (A  discussion  of  the  concept  of  ion  mobility  may  be  found  in  Ref¬ 
erence  1.)  Then,  since  the  total  velocity  of  the  charge  (ions)  is  equal  to  the  sum  of  the 
gas  motion  and  the  ion  motion,  the  following  relationship  exists: 

.  v^  ■  velocity  of  ions  relative  to  the  gas  +  velocity  of  the  gas. 

The  velocity  of  the  ions  relative  to  the  gas  is  equal,  to  KE.  Therefore, 

7,  =  KE  +  7 

where  v  is  gas. velocity.  The  current  density  becomes 

,  J 

=  Pg  I  KE  +  V  )  (4) 

using 

.  A  =  «  (V-'e  ) 


(1) 

(2) 

(3) 
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(5) 

Considering  the  relative  magnitudes  of  KE  and  v  for  a  typical  case  of  corona  wind: 

„  ~  2  cm^ 

“  volt-sec 

E  a  lOOQO  volts  _  5000  volts 

“  2  cm  spacing  ~  cm 

KE  =  100  m/sec 

The  usual  air  velocities  in  a  corona  discharge  are  below  5  m/sec.  Thus,  v  is  considerably 
smaller  than  KE.in  this  case  and  will  be  neglected.  If  air  is  blown  past  the  corona  point  at 
higher  velocities,  then  v  could  become  significant  (Ref.  6).  Therefore: 

J  a  p  Ke  s  K  ^  I  V  ■  7  >  7  (6) 

In  the  stagnation  region  under  the  impinging  corona  wind  stream,  velocities  will  be  very 
small.  The  boundary  layer  thickness' in  this  region  will  likewise  be  very  small,  and  the 
pressure  across  the  boundary  layer  will  be  approximately  constant.  Consequently,  the 
pressure  rise  due  to  the  electrical  field  Mil  occur  in  the  region  of  the  air  where  essen¬ 
tially  nonviscous  potential  flow  occurs.  Consequently,  in  the  Navier  Stokes  equations 
(Refs,  9,  10): 


and 

T  =  e  (^-7)  (  Kr+T  ) 


=  Xj  - 


X  =  F  ; 


,  M  dA 

-  -h  ii  V  ui  +  -  - 

axi  ^  '  3  a,. 

auj 

17“  ’ 


where  F,  the  elearical  body  force,  is  equal  to  X,  the  body  force  exerted  on  the  fluid.  The 
pressure  gradient  and  the  electric  field  body  forces  are  of  considerably  greater  magnitude 
than  the  viscous  or  inertia  forces,  and  consequently  the  Navier  Stokes  equations  reduce 
to: 


V  P  =  iP.  E 


and  . 


Vp  - 


follows.  This  indicates  that  pressure  gradient  Is  proportional  to  current  density  and  in¬ 
versely  proportional  to  ion  mobility.  In  a  one-dimensional  case.  Equation  (9)  becomes 

a  p  ^  j  • 
a  y  ’  K 
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^0 

In  a  one -dimensional  model  of  cross-section  A,  the  current  density  would  be  related  to 
the  current,  i,  by: 


Therefore; 


This  predicts  a  linear  rise  in  pressure  as  either  the  current  or  the  electrode  spacing  is 
increased,  is  a  reference  distance  located  a  very  short  distance  away  from  the  region 
of  intense  corona  discharge. 

If  the  electrode  geometry  is  cylindrical,  then  for  a  constant  current. 


and  A  is  given  as  a  function  of  radius  by: 

A  ?  ^  f 


where:. 

I  is  the  length  of  the  fine  wire  electrode, 

9c  is  the  angle  subtended  by  the  outer  electrode, 
and 

r  is  the  radius.  . 

Therefore: 


kX  0c  ''’’o 


dr 


f  '  dr 


(12) 


If  the  outer  electrode  radius  is  R,  and  the  reference  . radius  at  the  edge  of  the  intense 
corona  is  Tg ,  then  the  pressure  at  R  is; 


(13) 


This  relationship  would  indicate  that  the  pressure  is  once  again  proportional  to  the. current, 
but  varies  as  the  logarithm  of  the  radii  ratio  as  either  the  spacing  or  the  wire  size  is 
changed. 
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If  relationships  between  field  strength  and  current  are  desired,  then  it  can  be  shown 
(Refs.  2,  8)  for  the  linear  case: 


2i 

e  KA 


( y-  Vo* 


where; 


(14) 


^0 


For  the  cylindrical  case,  the  relaticMishlp  can  be  shown  to  be: 


(  — )  + 
'  «K  ' 


2 


.“r  •  '  '0  *“0 

where  E  is  the  field  strength  at  r  (Ref.  2). 


(r„  -21  ^ 


(15) 

(16) 


A  complete  listing  of  the  various  relationships  of  current,  voltage, and  pressure  for  simple 
geometries  is  given  by  Steutzer  in  Reference  7. 

Cobine  (Ref.  2)  derives  the  relationship  between  corona  current  and  applied  voltage  for 
a  cylindrical  electrode  configuration  in  the  form 


2KV  (V-Vc) 
R*  i  n 


(17) 


where  Vj  is  the  voltage  at  which  the  corona  starts. 

He  points  out  that  the  general  relationship  between  current  and  voltage  for  corona  dis¬ 
charge  is  of  the  form 


i  -  CV  (  V  -  Vg  )  (18) 

where  C  is  a  ccxistant  which  depends  primarily  uprai  the  geometry  of  the  system.  This 
form  has  been  checked  by  many  investigators  and  has  been  found  to  hold  quite  well  at 
low  current  values.. 

To  determine  the  pressure  rise  between  the  fine  wire  and  flat  plate  configuration  used 
in  the  experiments.attempts  were  made  to  derive  the  relationships  between  the  space- 
charge-limited  current,  voltage,  and  pressure.  The  boundary  conditions  in  the  problem  of 
a  round  cylinder  at  one  side  and  a  flat  surface  at  the  other  side  proved  difficult  to  handle 
mathematically.  Since  it  was  desirable  to. have  explicit  equations  for  the  field  and  current 
distribution,  numerical  solutions  and  representations  were  not  attempted.  Approximations 
based  upon  a  cylindrical  configuration  were  tried,  but  the  correlation  between  these  ap¬ 
proximations  and  the  pressure  curves  from  the  test  data  was  very  poor. 
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Pressure  Rise  with  Low  Charge  Density 


Since  small  currents  had  produced  very  significant  results  in  the  experiments, the 
possible  use  of  the  space-charge-free  field  equations  was  investigated.  Without  space 
charge,  the  problem  of  determining  the  field  becomes  elementary.  The  case  of  a  long 
wire  above  an  infinite  conducting  plane  is  covered  in  many  sources  (Refs.  11,  12,  13). 
For  close  spacing  and  small  wire  size,  it  was  assumed  that  the  test  configuration  might 
be  represented  adequately  by  the  following  theoretical  model. 


The  ^uaticsis  for  the  two-dimensional  field  based  upon  tbe  method  of  images  as  adaptec^ 
from  Reference  11  are: 


2  ir€ 

X 


2  IT 


r  ! 


y  +  o  ■ 


y— o 


(19) 


(20) 


where  the  various  distances  are  represented  in  Figure  36  and  x  represents  charge  per 
unit  length. 


At  the  surface  of  the  plate  y  =  0,  and  R  ^  ■  R^  =  R.  This  leads  to  the  following  relations: 
E,  =  0  at  the  plane  surface 


and 


X 

irt 


where  R®  *=  a®  +  x®. 

The  equation  for  voltage  is: 


V 


X 


which  is  valid  for  large  spacings  relative 


to  wire  size.  Therefore: 


(21) 


(22) 


Letting  Z  =  x/a 


2  (t)  r  1 

X*  +  0®  -I 

'  fa  ' 


2(1) 
in  (4^). 


+  Z 


(23) 


(24) 


The  relationship  for  the  pressure  in  a  dielectric  due  to  the  applied  electrostatic  field 
(Refs.  14,  15)  is  given  by: 
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A  p  =  *  E*  (25) 

Strictly  speaking,  this  relationship  holds  only  for  the  pressure  rise  in  a  dielectric  between 
two  parallel  plates.  It  does  not  hold  precisely  for  a  corona  discharge.  However,  Stuetzer 
{Ref.  7)  indicates  that  relationships  of  the  form 

Ap  =  -(G) 

hold  quite  well  for  corona  discharges.  G  is  a  constant  which  depends  on  geometry.  The 
maximum  value  for  E  will  occur  near  the  plate  where  E,  is  small.  Consequently,  as  an 
approximation 

Ap  =  «  E®  (6)  (26) 


with  the  constant  (G)  left  undetermined  for  the  moment. 

Using  this  relationship  and  Equation  (25),  the  expression  for  pressure  can  be  written; 


When  X  ■  0,  z  »  0,  pressure  at  the  centerline  can  be  expressed  by: 


where  p^  is  centerline  pressure. 

If  for  convenience,  the  reference  atmospheric  pressure  is  assumed  to  be  zero,  then 
APc  *  Pc .  Therefore: 


P 


(29) 


This  expression  provides  a  pressure  distribution  over  the  flat  plate,  assuming  there  are 
no  space  charge  effects,  that  p  =  I  c  E®  (G)^  and  that  the  influence  of  Ex  is  relatively  small. 

Equation'(29)is  plotted  in  Figure  37.  To  check  whether  this  expression  would  provide  a 
reasonable  representation  of  the  actual  phenomena,  it  was  compared  with  the  experimental 
data  and  the  results  are  plotted  in  Figures  38,  39,  and  40.  These  figures  present  the  same 
pressure  distribution  data  as  was  presented  in  Figures  30,  31,  and  32,  except  that  the  data 
have  been  normalized  by  dividing  the  magnitudes  of  all  pressures  by  the  centerline  peak 
pressure.  The  distances  along  the  plate  were  transformed  by  introducing'a  new  variable, 

Z  (Z  =  x/a,  and  a  is  the  distance  from  the  wire  to  the  plate).  These  plots  show  that  thfe 
form  of  the  pressure  distribution  for  a  given  wire  spacing  remains  approximately  the 
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same  as  the  current  is  varied  widely.  In  addition,  the  shape  of  the  distribution  remains 
essentially  the  same  with  changes  in  the  spacing.  The  normalized  distribution  is  more 
peaked  with  the  6-cm  spacing,  but  the  differences  are  not  large.  Large  changes,  however, 
can  arise  from  small  zero  shifts  in  the  nriicromanometer.  These  effects  are  shown  in 
Figure  41^which  presents  the  normalized  data  of  Figure  33.  Data  obtained  at  the  lower 
three  values  of  current  had  small  random  shifts  in  the  zero  reference  of  the  manometer. 
As  can  be  seen  from  Figure  41,the  wide  variance  makes  these  data  unusable. 


As  can  he  noted  from  Figures  38,  39.  and  40,  correlation  between  the  shape  of  the 
pressure  curves  for  the  test  data  and  for  theory  is  quite  good  at  smaller  values  of  a 
at  larger  values  of  a.  however,  correlation  becomes  poorer.  This  may  possibly  be 
attributed  to  the  fact  that  the  assumed  infinite  plane  of  the  two-dimensional  model  be¬ 
comes  a  poor  approximation  as  the  spacing  becomes  large.  Based  upon  this  comparison 
with  the  test  data,  the  variation 


is  considered  a  reasonable  representation  so  long  as  a  is  considered  small.  The  vari¬ 
ation  provides  information  on  the  shape  of  the  pressure  distribution  only. 

In  attempting  to  calculate  pj  from  the  form 


(6  )  , 


it  was  found  to  be  difficult  to  determine  the  geometric  constant  G  because  of  the  boundary 
conditions  mentioned  in  the  space-charge  case.  Consequently,  the  data  obtained  during  the 
pressure  tests  were  used  in  calculating  center  pressures.  The  relationship  between  center 
line  pressures  and  current  was  in  accord  with  the  equations  for  the  one-dimensional  case 
and  the  cylindrical  configuration. 


One-dimensional:  p  =  — i —  (  y  _ 

Cylindrical;  p  *  _L —  Xn  {  —  )  (13) 

K  ^  ' 

The  test  data  shown  in  Figures  27,  28,  and  29  verify  the  linear  relationship  of  pressure 
and  current.  Using  p  =  c  i  i,  the  constant  Cj  becomes  essentially  a  geometric  shape  factor 
for  a  given  ion  mobility.  Ci  then  is  determined  from  Figures  27.  28.  or  29  fdr  the  specific 
wire  size  and  wire  spacing  chosen. 

Thus,  in  establishing  the  mathematical  model  of  the  pressure  produced,  an  experimental 
constant  has  been  introduced.  The  resulting  expression,  however,  has  the  advantage  of 
being  quite  simple  in  form. 

CONVECTIVE  HEAT  TRANSFER  WITHOUT  APPLIED  ELECTRIC  FIELDS 

A  brief  summary  of  relevant  heat  transfer  problems  will  be  given  to  provide  a  back¬ 
ground  for  this  work.  Flat-plate  forced  convection,  vertical-plate  free  convection,  and 
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stagnation-point  heat  transfer  will  be  covered.  The  starting  point  for  solving  all  three 
cases  is  with  the  following  equations  (Ref.  9); 

Continuity: 


d  U| 
d*i 


*  0 


Navier  Stokes: 


Energy: 


duj 
d  r 


=  Xi  - 


—  +  /iV  uj  +  - - — 

dxj  3  dxi 


(30) 


(7) 


d  Ca  T  I  Hn  3  3t  M  I 


dr 


Jm  dr  d*j  dxj 


■»(« 


(31) 


where: 


(7f)‘ 

t  (  is.  .  ±L.)‘  .  (±i 

'  3y  .  dz  '  '  dx 

q^  “  radiation  energy 

q^  =  internal  heat  generated 

Xj  =  a  body  force 

r  =•  time 

dui 


d  w 

/  du 

d  z 

)  * 

V  dy 

dw 

2 

-  A®  1 

dx 

-  ■ 

) 

3 

•  A  ) 

A  = 


d  Xi 


(32) 


Forced  Convection-Laminar  Fiow 


With  the  usual  assumptions  of  Biasius  flow  (Refs.  4.  16).  the  equations  for  the  case  of 
heat  transfer  to  a  flat  piate  in  iaminar  flow  become 

I 

.  Continuity: 


Navier  Stokes: 


20  . 


d  u  d  V 

d  X  d  y 


0 


d  u 

u  - 

dx 


+  V 


du 

dy 


V 


d\ 
d  y* 


(33)- 


(34) 
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Energy; 


dT  dT 

u  -  +  V  - 

d  »  dy 


Use  of  the  stream  function 


and  the  transformation 


'J'  =  (  1/  X  U  )*  f  (7)) 


7)  =  y  {-^^Y 

•allows  the  Navier  Stokes  equations  to  be  transformed  to  an  ordinary  differential  equation 

•  • 

f  (Tf)  riTj)  =0  (37) 

the  solution  for  which  is  tabulated  in  many  sources  (Refs.  4,  10).  The  energy  equation  is 
handled  through  the  same  transformation  and  becomes: 


^  iv  -  (■ 


'^m  Cp 


•)  Npr 


The  solution  to  this  equation  is  given  by; 


■-'{ir-f 


A  reasonable  approximation  to  this  solution  (Ref.  17)  is  given  by: 


=  k 


where  Np^  is  the  Prandtl  Number. 

Free  Convection  Vertical  Plate 

In  the  free  convection  from  a  vertical  plate,  a  body  force  arises  due  to  the  effect  of 
gravity.  The  equations  in  this  case  become  (Ref.  9); 

Navier  Stokes; 


Energy: 


P  (  u  +  V  )  =  X  -  +  u 

dy  d  X  d  y^ 


dT  dT  k  d  T 

u  -  +  V  -  =  -  - 

dx  d  y  p  Cp  dy‘ 


The  transformation  of  the  equations  in  this  case  uses 


ASD-TDR-62-650 


where 


'I' 


cy 

TT 


s_ 

4 


Fiiy) 


^o  “  stream  temperature 

=  wall  temperature 

9  .=  T-T  /T  -T  . 

■  o'  w  o 

The  equations  then  become: 

Navier  Stokes; 


+  3  FF"  -  F‘  +0  =  0 


Energy;  . 


(43) 


9"  +  3  Np,  FS’ s  0 

Solutions  for  this  set  of  equations  are  listed  in  References  4  and  10. 
Stagnation -Point  Heat  Transfer 


(44) 


The  stagnation-point  problem- arises  when  a  fluid  stream  impinges  upon  a  surface  and  a 
zone  of  zero  velocity  is  formed  under  the  stream.  The  case  of  an^finfr^stream  striSn^ 
a  cylinder  or  flat  plate  is  well  covered  in  the  literature  fRefs  0  Rntv.  .vi  • 
experimental  evidence  support  the  assertion^sSTthe 

variaS^toli  onlv^nfhe'^-^®'''^  stagnation  point.  This 

variation  holds  only  in  the  vicinity  of  the  stagnation  point,  but  it  is  most  useful. 

equaS tcSe?'*''"  ^  ^^rm  in  the  Navier  Stokes 


-  - ^  -  i3 

O  i  ^1 


(45) 


and  the  Navier  Stokes  equations  become: 


d  u 


d\» 


TT'"''  ~ 


dy 


-£JL 

d  y* 


(46) 
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Using 

j, 

'1'  =  (  V  /8,  )*  X  FCij) 

.  /  \T 

the  equation  is  transformed  to  the  ordinary  equation 

F*"  -I-  FF"  _(F*)*+  1=0. 

The  energy  equation  for  this  case  is  given  by 

dT  .  dT  k  d*T 

u -  +  w  -  =  -  - - 

^x  dy  /)Cp 

Which,  with  the  transformation,  becomes 

6“  -H  Np^  F  e*  =  P  (49) 

Where 


Td~ 


The  solution  of  these  equations  is  given  in  Reference  16. 

The  foregoing  brief  sketch  of  three  of  the  fundamental  problems  of  convective  heat 
transfer  are  included  solely  as  background  material.  Some  of  the  efforts  made  in  tackling 
the  problem  of  heat  transfer  with  corona  discharge  were  based  upon  the  concepts  involved 
in  these  solutions.  Their  Inclusion  provides  for  ready  reference. 

COMBINED  CORONA  WIND-CONVECTION  PROBLEM 

In  order  to  formulate  a  suitable  mathematical  model,  it  is  necessary  to  set  down  as 
completely  as  possible  the  physical  system  to  be  described.  In  the  present  work,  two 
separate  physical  cases  will  be  considered.  The  corona  wind  impinging  on  a  horizontal 
plate  with  the  bottom  heated  constitutes  one  case,  and  the  corona  wind  impinging  on  a 
heated  vertical  plate  constitutes  the  second.  These  two  systems  are  shown  in  Figures 
42  and  43.  In  the  first  case,  the  bouyancy  force  is  symmetrical  about  the  center  of  the 
plate.  Any  streaming  due  to  thennal  bouyancy  or  corona  wind  is  additive  and  no  regions 
of  stagnation  should  exist. 

The  bouyancy  force  in  the  case  of  the  horizontal  plate  acts  somewhat  differently  from 
that  for  the  vertical  plate.  For  a  vertical  plate,  the  bouyancy  force  acts  along  the  plate, 
and  causes  the  convective  streaming  to  move  along  the  surface.  For  a  horizontal  plate, 
the  bouyancy  force  acts  perpendicular  to  the  plate,  and  any  streaming  occurs  because  of  ■ 
the  upward  movement  of  the  heated  air  at  the  edge  of  the  plate.  Consequently,  in  the  center 
of  the  plate,  the  bouyancy  force  which  acts  perpendicular  to  the  small  streaming  taking 
place  will  be  neglected  in  the  analysis  of  the  horizontal  plate. 


(47) 

(48) 
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The  vertical  plate  presents  an  interesting  but  formidable  problem.  The  same  complex 
additive  corona  effects  and  thermal  effects  occur  as  in  the  first  case,  and  in  addition,  the 
magnitude  of  corona  wind  will  approach  that  of  the  free-convection  streaming  at  certain 
values  of  current.  Thus,  at  least  two  separate  problems  present  themselves.  In  the  region 
above  the  centerline  of  the  impinging  jet  of  air,  the  two  streams  will  be  additive  and  a 
combined  free-convection  and  forced-convection  problem  results.  The  thermal  body  force 
in  this  case  is  in  the  same  direction  as  the  external  corona  wind  stream.  In  the  region 
below  the  centerline, the  situation  is  changed.  Here  the  corona  wind  streams  in  the  opposite 
direction  to  the  thermal  body  force.  At  some  value  of  current,  the  magnitude  of  the  corona 
wind  will  be  approximately  the  same  as  that  of  the  upward  streaming  heated  air.  Although 
it  is  doubtful  that'  the  distribution  of  velocity  in  the  two  streams  as  they  impinge  will  be 
the  same,  the  net  result  is  that  some  type  of  stagnation  region  should  result.  This  com¬ 
plex  problem  is  illustrated  graphically  by  the  interferometer  pictures  of  the  vertical 
heated  plate  shown  in  Figure  13. 

The  approaches  used  in  attempting  to  solve  these  problems  included  consideration  of 
solutions  of  the  Navier  Stokes  and  energy  equations  with  the  concept  of  an  impinging  two- 
dimensional  finite  jet,  and  the  use  of  the  Von  Karman  integral  equations. 

Two-Dimensional  Jet  Approach 

Study  of  all  the  data  and  photographs  of  the  corona  wind  indicated  that  the  cross  section 
of  the  impinging  jet  of  air  appeared  to  be  very  narrow  before  it  was  deflected  by  the  plate. 
Pictures  of  both  the  horizontal  plate  and  the  vertical  plate  at  high  currents  indicated  that 
the  streaming  after  impingement  was  approximately  symmetrical.  A  possible  model  then 
was  conceived  to  be  a  thin  two-dimensional  jet  of  air  impinging  upon  a  flat  surface.  The 
case  of  the  horizontal  plate  with  the  thermal  body  force  included  was  considered  from  this 
viewpoint  first.  If  this  approach  would  turn  out  to  be  promising  it  was  felt  that  the  width 
of  the  jet  could  be  determined  from  the  test  data. 

The  case  of  the  horizontal  flat  plate  with  an  impinging  jet  was  considered  in  the  follow¬ 
ing  manner.  As  a  first  approximation,  it  was  assumed  that  the  viscosity  and  boundary  lay¬ 
er  would  not  markedly  influence  the  flow  field  and  pressure  distribution  of  the  jet  of  air. 
This  assumption  allowed  the  use  of  potential  flow  theory  in  the  region  away  from  the  sur¬ 
face  of  the  plate.  One  possible  technique  for  handling  this  case  is  by  using  the  hodograph 
transformation  (Refs.  18,  19).  The  problem  is  shown  in  Figure  44. 

Using  the  transformation 

f  .  (u-iv)  /rm 


the  jet  is  mapped  onto  the  upper  half  of  the  unit  circle. 

maps  the  upper  half  of  the  unit  circle  onto  the  upper  half  of  the  t  plane.  The  flow  in  the  t 
plane  is  equivalent  to  a  source  flow  of  strength  . 

The  complex  potential  w(z)  =  u  +  i  v  becomes,  for  this  case: 
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=  Jil-  in  (-^) 

2ir  V  t®-l  / 

Jjti 

C*-i 


=  4fLi. 

2.ir  ,  '  r*_  I  ' 


dw  dw  dC 

Using  3  — ^  - ,  it  can  be  shown  that 

dz  df  dz 


dz 


=  ii-  /_J _ ! _ \ 


from  which 


—  =  (-^  -  )  (ton  't  +  tonh't  ). 

b  w  . 


At  the  surface  of  the  plate,  v  =  0  and  this  equation  becomes 


JL  (4-)  JL 

b  \  ir  /  \  U 


4“  ton  h 


- 1 


u  / 


(51) 


(52) 


(53) 


This  is  an  implicit  equation.  Expressing  the  terms  in  series  and  reducing  the  equation. 


Using  a  reversion  of  series  technique 


The  form  of  the  equations  for  velocity  and  pressure  calculated  on  an  incompressible  basis 
is  shown  in  Figure  45. 

The  most  revealing  aspect  of  this  figure  is  that  the  linear  term  of  the  expression  for  u/.  j 
holds  quite  well  in  the  vicinity  of  the  centerline.  If  the  pressure  and  velocity  relationship 
of  the  corona-wind  jet  would  be  similar  to  these  curves,  then,  potentially,  a  very  simple 
calculation  for  boundary- layer  profile  and  heat  transfer  could  be  made.  In  the  central  re¬ 
gion  where  x  <  —  b  ,  the  stagnation  solution  discussed  previously  would  apply. 

4u 

and  in  the  outer  region  where  x  >  -  b  .  the  flat-plate  laminar  flow  solution  would 

apply. 


A  comparison  between  the  theoretical  thin- jet  pressure  distribution  and  a  typical  pres¬ 
sure  distribution  obtained  from  test  is  shown  in  Figure  46.  Several  attempts  were  made 
to  adjust  the  relative  shapes  of  the  curve  by  choosing  various  values  of  b,  the  theoretical 
jet  half-width.  Since  the  form  of  the  pressure  curves  did  not  agree  well  with  the  test  curves, 
it  could  be  expected  that  correlation  with  heat  transfer  coefficients  based  upon  this  method 
would  also  not  be  very  good. 
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Navier  Stokes  Approach 

Since  both  the  test  data  and  the  thin-jet  theory  indicated  that  large  pressure  variations 
would  be  confined  to  a  region  close  to  the  centerline,  it  was  felt  that  the  combined  thermal 
body  force-corona  wind  problem  could  be  tackled  in  the  outer  regions  of  the  plate.assuming 
the  wind  velocity  in  this  region  to  be  approximately  constant.  The  physical  model  which 
was  considered  next  was  that  of  the  second  case  as  represented  by  Figure  43.  This  figure 
illustrates  the  two  regions  of  flow,  above  the  centerline  and  below  the  centerline.  If  the 
free-convection  effects  are  large,  and  the  upper  half  of  the  plate  is  considered,  then  the 
body  force  acts  in  the  direction  of  external  flow.  In  this  case,  the  body  force  is: 


where  the  term  ( — - — )  ( — ^  is  dropped  from  the  energy  equation  because  the  velocities 

'JmCp  '  '  d  y  ' 

are  extremely  small.  This  set  of  equations  is,  of  course,  exactly  the  same  set  as  was  used 
for  the  free-convection  vertical  plate.  The  significant  change  in  the  problem  comes  with 
the  different  boundary  conditions.  As  the  coraia  wind  streams  over  the  surface,  the  vel¬ 
ocity  distribution  will  vary  from  zero  at  the  surface  to  some  maximum  away  from  the 
surface  and  then  decay  back  to  zero  at  a  distance  farther  out  from  the  plate.  The  shape  of 
this  total  distribution  is  complex.  The  simpler  problem  of  the  spreading  of  a  free  jet  is 
even  difficult  to  handle.  In  the  case  of  corona  discharge,  it  is  difficult  to  know  the  form 
of  the  highly  localized  velocity  distribution,  and  consequently  the  shape  of  the  wake  is  in 
doubt.  This  shape  would  form  the  outer  region  of  the  corona  wind  after  it  turned  the  comer. 
On  the  other  hand,  the  width  of  the  boundary  layer  growing  beneath  the  corona  wind  would 
be  very  small  for  some  distance  away  from  the  center  stagnation  region.  It  was  assumed, 
consequently,  that  the  thickness  of  the  corona  wind  would  be  large  relative  to  the  boundary 
layer  thickness,  and  that  the  velocity  at  the  edge  of  the  boundary  layer  would  be  the  velocity 
of  the  peak  of  the  corona  wind.  These  assumptions  are  represented  in  Figure  47  where 
schematic  representations  of  temperature,  and  velocity  profiles  are  shown.  ,  the  imping¬ 
ing  velocity,  is  the  velocity  of  the  corona  discharge. 


Using  the  free-convection  transformations,  equations  (56)  and  (57)  become: 


F'"  +  3  FF" -2  (F'l*  4-  0  =  O  (58) 

0"  4-  3  Np,  F  0'  =  0  (59) 
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The  boundary  conditions  are; 


which  lead  to: 


y  =  0,  If  =  0  ,  0=0,  V 
y  ^8,  17-oc,  0=  ug,  V 


0  ,  T  =  T„ 
0.  T=  To 


F=0,  F'  =  0,  fi=l,  0  y=0 


F"=0,  F' 


4  1/  c®  K-i- 


®  y  —  S 


These  relationships  imply  that  the  boundary  conditions  change  at  every  position  up  the 
wall.  If  this  approach  is  used,  solutions  for  the  equations  would  have  to  be  made  in  dis¬ 
crete  steps  for  selected  values  of  x.  , 


For  the  case  with  the  body  force  acting  in  the  direction  opposite  to  the  direction  of  flow 
(Ref.  20),  the  sign  of  the  body  force  term  is  changed  and  the  equations  become 

F“'  +  3FF"  -2(F‘)®  -0  =0  (60) 

e"  +  3Np^  FB'  =  0  (61) 


The  boundary  conditions  are  the  same  as  for  the  previous  case.  Consequently, a  solution 
proceeding  in  this  fashion  would  have  to  be  done  in  discrete  steps  also. 

Different  forms  for  the  transformations  were  tried,  but  no  clear-cut  approach  for  handling 
the  combined  free  convection-forced  convection  was  evident.  Since  a  review  of  other  com¬ 
bined  convection  work  did  not  reveal  anything  directly  applicable,  it  was  decided  to  follow 
the  integral  equation  approach. 

Von  Karman  Integral  Equation  Approach 


The  Integral  equations  of  boundary  layer  and  heat  transfer  theory  for  steady  flow  (Ref, 
10)  are  given  as: 

Momentum: 


-g 


H- 


0 


(62) 


Energy: 


fn  p 


where  6^  is  the  thickness  of  the  thermal  boundary  layer.  Because  the  velocities  involved 

in  the  phenomenon  are  so  small,  the  viscous  dissipation  term  was  dropped  in  the  energy 
equation. 
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The  procedure  used  in  solving  the  integral  equations  is  to  assume  reasonable  forms 
for  the  velocity  and  temperature  distributions  in  the  boundary  layer  which  meet  certain 
prescribed  boundary  conditions.  Although  the  technique  is  not  precise,  it  does  offer  insight 
into  the  problem  as  well  as  provide  solutions  useful  in  engineering. 

ANALYTICAL  SOLUTION  OF  THE  HORIZONTAL  PLATE  CASE 

In  the  application  of  the  Von  Karman  integral  equations  to  the  case  of  heat  transfer 
from  the  bottom  of  a  heated  flat  plate  subject  to  corona  wind,  it  was  assumed  that  the 
thermal  body  force  was  small  and  could  be  neglected.  The  integral  equations  then  become: 


,8r 


,^87 


St 


wo,  -|t (^)  ■ 


dT 

dy 


(65) 


In  this  case,  it  was  decided  to  attempt  to  determine  what  variation  in  boundary-layer 
thickness  and  heat  transfer  coefficient  could  be  obtained  throughout  the  entire  range,  from 
the  centerline  outwardly.  The  body  force  effects  due  to  the  applied  electrical  field  were 
included  through  the  use  of  the  pressure  rise  distribution  previously  derived  under  the 
section  on  electric  field-corona  discharge  relations.  Thus, the  analysis  consisted  pri¬ 
marily  of  determining  the  variations  of  boundary-layer  thickness  and  heat  transfer  co¬ 
efficients  under  the  influence  of  a  prescribed  pressure  distribution. 


Pressure-Velocity  Relation 

In  considering  the  fluid  flow  in  the  regicai  of  the  impinging  jet  of  corona  wind,  it  is 
necessary  to  examine  the  nature  of  the  relationship  between  velocity  and  pressure.  For 
incompressible  flow,  Bernoulli’s  equation  holds.  If,  however,  a  body  force  is  introduced, 
the  relationship  Is  not  as  simple,  even  if  the  viscous  effects  are  considered  small.  From 
thermodynamic  reasoning  it  can  be  shown  that  if  work  is  done  on  the  fluid,  for  example 
by  the  body  force,  then  the  constant  of  integration  will  change  across  streamlines,  and 
the  condition  of  irrotationality  will  not  hold  even  for  incompressible  flow.  Because  of 
this,  one  cannot  assume  the  relationship  Pc  =  P  +|  pu®  will  hold  throughout  the  fluid 
stream  above  the  boundary  layer  without  more  careful  consideration. 


In  the  region  outside  the  boundary  layer,  the  viscous  effects  can  be  considered  small 

(66) 

(67) 


and,  consequently,  the  Navier  Stokes  equations  become 

0  ( 

du  du  1 

i  -  X  - 

P  \ 

U  '  ^  V- 

dx  dy  ‘ 

^  ax 

0  ( 

dw  dv  \ 

1  =  Y  - 

d*  dy  * 

'  ay 

If  ■  Pc  •  using 

dxj 

where  V  is  the  potential 

(voltage),  then 
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This  condition  holds  for  all  values  of  u  and  vif  7®u  =  0,  V®v»0,  and  implies  that 
if  pc  can  be  considered  a  scalar  constant,  then  irrotationality  will  follow.  The  foregoing 
demonstrates  that  the  use  of  the  condition  of  irrotationality  is  reasonable  for  the  case  of 
very  small  space  charge.  Therefore; 


( u  w  \  dK  =  -  d  (  p  +/5  V  ) 

V  ax  d%  f  p  c 

/j_  (  u*  -(-  V*  )  -  -  P  -  Pc  ^  -i-  Po 


(76) 


p  +  p_(u*+v*)  =  p— /pV  (77) 

where  po  is  the  constant  of  integration.  Near  the  flat  plate  surface,  v  =  0,  and  V  approaches 
zero.  Consequently,  Bernoulli’s  equation  holds  approximately  near  the  wall: 


P  +  P—  u*  =  p  .  (78) 

2  0 

When  X  =  0,  u  =  0,  p  =  Pc  ,  it  follows  that  po  =  Pc  ,  where  pc  is  the  pressure  at  the  center- 
line  of  the  corona  discharge. 


Using  the  reiation  previously  derived  for  the  pressure  variation  under  a  corona  discharge. 


the  Bernoulli  equation  becomes; 


f— L_]“ 


(29) 


(79) 


With  this  relation,  the  soiution  of  the  integral  equations  can  be  undertaken.  Three  forms 
of  velocity  and  temperature  profiies  are  assumed:  linear,  parabolic,  and  cubic.  Because 
the  parabolic  distribution  soiution  provides  the  clearest  picture  of  the  techniques  involved, 
it  will  be  carried  through  completely.  The  linear  and  cubic  solutions  will  be  indicated 
oniy. 


Paraboiic  Soiution 


The  parabolic  profiie  assumed  is; 


which  meets  the  boundary  conditions  of: 

y  =  0  ,  u  =  O 

y  =  8  ,  V  =  ux 


du 

dy 


=  0 


S  m  T  -Tw 


e  = 


(80) 


(81) 
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which  meets  the  boundary  conditions  of: 

y=0.  T=T^.  8  =  0 

y  =  St  T  =  To  .  9=9^ 


Substitution  of  the  velocity  profile  into  the  momentum  equation 


°  u  d  ^  s  dp  t  <>»  \ 

J  ^  ,  J  ^ody  =-S  —  -  ^  (1T)w  = 


du 


dx 


dx 


d  y  'y  =  0 


leads  to 


15 


(u„S)  -t-S^ 
*  dx 


(64) 


(82) 


which,  upon  further  manipulation,  becomes 


dS*  r  3  r  *  1  ^ 

'] 

30  1/ 

(83) 

dx  L  0^2  dx  '  *  ■  / 

“x 

Considering  the  bracketed  term  and  the  Bernoulli  relation 

z 

Uk 

p  +/>  —  ■  p. 

(84) 

X  p  p 

(85) 

_i _ 1_  /  *  ^  _5p_  \  _  3  d 

Ux*  dx  '  *  pi  Uj*  dx 

(86) 

Since  pc  is  a  constant  with  respect  to  x 

3  d  /  2Pc  3p  \  _  9 

1 

dp 

(87) 

Uj*  dx  '  />  P  '  2 

//>Ux*' 
'  2  ‘ 

)  ■ 

But 


_dP _ _  /  P'>n  \ 

dx  '  dx  '  2  ' 


(88) 


and 


9 

1  ^  \  2  / 

2 

^Ux*  dx 

2  dx^  ^  2  ^ 

(89) 
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The  momentum  equation  becomes: 


d8* 

+  -i- 

r  J-  1 

8*=  . 

(90) 

dx 

+  ■ 

2 

^  dx 

k  2 

The  solution  to  this  equation  is 


For  the  boundary  layer  thickness  to  remain  finite  at  x  =  0.  u,  =  0.  c  must  be  zero. 


(93) 


(79) 


where  u  =  u* ,  the  velocity  outside  the  boundary  layer. 


8  ^  =  15  /i  a 


where 


(94) 


z  =  x/a 

a  =  spacing  of  corona  wire  from  plate. 

Using  a  somewhat  tedious  process,  the  integral  for  this  parabolic  case  can  be  integrated 
exactly,  y*  is  expanded  and  integrated  term  by  term  to  give: 
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(2)  =  /  / 


dz  =  2  —  0.  899023 


■+  ton  2 


+  0.7333984  - O.  58 6 71 88  - - ^ 

(1+2*)  (1+2  ) 


0.35424  tl 


(H-2*) 


2^4 


-  O.  314  8810 


(1+2*) 


2  18 


+  0.0  7738  10 


(1  +  2*) 


2x8 


+  0,0714286 


(1+2*) 


2l7 


(95) 


Thus  an  explicit  solution  for  the  square  of  the  thickness  of  the  boundary  layer  is  obtained. 
The  behavior  of  the  expression  for  6®  in  the  neighborhood  of  the  origin  must  be  investi¬ 
gated  closely,  because  the  expression  y®/®  approaches  zero  at  the  origin.  Applying 
L’Hopital’s  rule,  let 


f  (2) 
g  (z ) 


/y-*d2 


y  * 


Then 


f'(2l 
g'(2  ) 


f 


2 

9 


iz 

dz 


_2 

9' 


4  2(1+ 2®  r’  '® 


'  (2  +  r*)*  (1+2*)* 


(96) 


lim  - 
2—0 

At  the  origin  then 


f  (2) 
g  (z) 


18 


8o 

For  larger  values  of  z 


S*  =  15/4  0  ^  ^  ■ 
The  shape  of  the  boundary  layer  is  given  by 


(97) 


(98) 

(99) 


(100) 
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_  ' 

(15^.7^)+  [  r‘  . 


(101) 


This  function  is  plotted  in  Figure  48.  The  function  shows  the  finite  value  at  the  origin,  and 
the  subsequent  increase  as  z  gets  larger. 

The  energy  equation  will  be  considered  next. 

ir{  (“ar)y=o' 

Substituting 
0  =  T  -  T* 


38 


dp 


^  [  r'(u0-ue„)  d 1 - ^ 

d  X  '■  *^0  ®  8t  P  d* 


N/ 


u  d 


St 


.„(^)  (102) 

'  dy  'yso 


-0- 


il  f 

2  a  do 

(103) 

6  >■ 

5  J 

'  dx  *  T  ^j^Cp  1 

^  3  -1 

St 

where  8* 

38t 
3  X 


and  is  assumed  ncx  to  vary  appreciably  with  x. 
dp 

I  du,*  _ 8«-3 


dx 


-  20 


.-3  dx  I  g* 

,-5  0JmCp®o  ^  '  8^(8r-5  )u 


(104) 


which  solution  is: 


=  e 


-f  f,  (x)  dx  Jf,  (x)  dx 
/* 


(x)  dx  c  e 


(x)dx 


where 


(105) 


f  I  (x)  dx  =  { 


dp 

8<t-3  ax 


r  I  du. 

=  J  — —  —  -20 - 

dx  8*-5  P'^m'^p^o 


fj  (x)dx  ® 


120  a 


8*(  8»— 5 )  u- 
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Therefore: 


/f,  (X) 


r  -Q  1 

/S*-3  X 

1 

r  dp 

/  .  ■■  u »  1 

'  dx 

xe 

V  S«-5  i 

P  ‘*iti^p  ®o 

^  dx 

/  8»-3  \  20 p 

^  S»-5  >  ^J,nCpSo 


r  20P  . 

/  e  f,  (x)dx  =  /u  '  S»-5  '  P  Jn»Cpfio  (-|) - 120_a - 

^  J  X  S*(S*-5)u, 


_  /  0*-3  \  20p 

— !i2_5! /  U_e  '  8x-5  '  /’>Jn»Cp0o  dx 

S«(8*-5)  •'  * 


Substituting 


l  +  z*  ■ 
8*  -  3 


-  (  C*-~3  \  20 

^  S*-5  '  PJmCpSo 


2  2  Pc 

“x  =  — ^  r 


and  noting  that  for  the  solution  to  remain  finite  at  the  origin,  c  =  0, 

2  -laoaoe  i.  Pc  (tTp) 

8t - /  go-;—  / 


8*  (8*-5)  V  -j-  r  0 

Using  the  first  two  terms  of  the  series  expansion  for  e“ : 

-  2  -  S 


e  =  I  +  X 


2/  3/ 


leads  to 


The  expression  for  6,^,®  becomes: 

^2  I20aa  [  '““i  Pc(  1  +  32  )  ]  i.  r 


^  2  120 go  L  '  “I  V  i  +  z! 

’  8*l8*-5)  yi^y 

P 


/%^[>  +  “.P,{-^)]  dz  (111) 
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This  is  a  complex  expression  for  which  no  exact  solution  was  found. 

The  expresslCHi  Pj  will  be  examined  hirther.  At  the  centerline  Pc  ■  /L  ^  »  where  U 
is  the  peak  impinging  velocity. 


where 


8»  -3  aopc 

8«  -5  P 


8*  -  3 
S*--  5 


^Ek 


lOU^ 

JmS  8 


JmCp  8^ 


S»  -3 
8»-5 


lONEk 


(112) 


To  determine  the  order  of  magnitude  of  aiPc .  consider  typical  values  of  the  various  terms. 


Typically: 


9o  =  20* F 

U  =  10  ft/ sac 

Cp  *  0.24  BTU  /  lb  'F 


“  8.32  X  10“* 

For  values  of  6*  from  0  to  3,  the  expression 
value  of  1/2  which  occurs  when  6,  =  1, 


6*  -  3 

- - =-  lies  between  3/5  and  0.  Assuming  a 

0,-5 


a,  Pg  =  4.  16  X  10 


Since  the  pressure  relation 


[^1' 


(29) 


varies  between  1  and  0  as  z  varies  from  0  to  infinity,  a,  p^  (  ^  ^  2  )  will  always  be 

smaller  than  4.16  x  10“®  for  6*  =  1.  Values  of  6*  above  2  are  unlikely.  Evaluation  of  6.J-/6 

at  specific  points  verified  this.  If  6*  becomes  greater  than  3  the  analysis  would  become 
questionable,  certainly  in  the  region  6*  =  5.  Therefore: 


a,  p  4.16  xio 
a  <  a 


<  ~  1.004 

for  all  values  of  z. 

For  the  small  pressures  involved  in  a  corona  discharge,  then. 
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o,  p 


s  I 


2  -120  a  a  .  ^ 

■  r.s.(S.-5i  VIS  •' 

/» 

The  value  of  6,j,®  as  z  -  0  is  indeterminate.  Applying  L'HopitaTs  rule 

IM.  .. 

g(il  Y 

±  ±  ' 
f '(r)  Y^  (2  +  1)  (l+i*) 

y'(i)  Jy  4z(l+z*)'®  4 

dz 


X-  im 


2—0 


f(z) 
g  (z ) 


At  the  origin 


4 


-30  a  0 
S  « (S«-  5  ) 


(113) 


(114) 


(115) 


(116) 


By  definition  6*  =  fij/g*  The  value  of  6*  at  the  origin  can  be  determined  as  follows  using 
the  values  of  6  and  6®  at  the  origin. 


szrfcr) 

(t) 


le  Npy~' 

8*  ( 8*  —  5) 


(117) 


For  air,  with  Np^  =  0.71,a  graphical  solution  yields  6*  =  2.048. 

For  values  of  z  >  0,  the  ratio  varies  in  a  more  complex  fashion.  For  small 

pressures,  this  ratio  becomes; 


8, 


-^**Pr 

S*(8*-5) 


/  y*  dz 
Y 

/y^dz 

2 
Y  * 


(118) 
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As  z  approaches  infinity,  the  bracketed  term  approaches 

1. 


and  the  ratio  approaches 


At  very  high  z. 


8 


t 

» 


Np;' 

8*  (8*  -5  ) 


(119) 


For  a  Prandtl  number  of  unity,  6*  =  1.  For  air,  6*  =  1.05.  Thus  a  variation  in  6*  does 
occur  with  a  change  in  z.  This  is  contrary,  of  course,  to  the  assumption  originally  made 
in  solving  the  integral  equations.  The  effect  of  the  change  in  6*  will  be  considered  in 
evaluating  the  over-all  results. 


The  heat  transfer  coefficients  can  be  calculated  from 


=  - 


k  dT 

^  y  y  =0 
2 


Using  0  =  T  -  T^  and  5  , 


2k 


Expressing  this  in  terms  of  the  velocity  boundary  layer. 


(120) 


(121) 


(2  P  Pc  I 

(30o/tH 

yf 

relationship  found  in  the  previous  section, 

Pc  *  * 

h  =  -li-  I _ 

*  8*  (30o  /i  )''2  r  F(z)  -  '4 

L  J 

Near  the  origin  6*  =  2.048 


2k 


8« 


Using  the  pressure-current 


At  the  origin 


2k 

2.048 


{a 


(122) 


(123) 


(124) 


where  Ci  is  an  experimental  and  geometric  constant.  This  equation  predicts  that  the  heat 
transfer  coefficient  at  any  station  should  vary  as  the  one-fourth  power  of  the  current. 


The  shapes  of  the  velocity  and  temperature  profiles  which  are  assumed  may  materially 
affect  the  magnitudes  of  the  boundary  layer  thickness  and  heat  transfer  coefficient  but 
should  have  relatively  small  influence  on  the  form  of  the  variation  in  heat  transfer 
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coefficient  with  distance.  In  order  to  get  an  estimate  of  the  changes  caused  by  profile 
shape,  two  other  profiles  were  assumed,  linear  and  cubic. 

Linear  Solution 


For  the  linear  solution  the  foliowing  curves  were  assumed: 


«-a.  —  ;  8  .  T-T. 


with  the  boundary  conditions  of 

y*0,  TsTj,  ,e=0,u=0 
y  s  S  ,  u  s  Ux 

»  =  St.  t  =  To  .  e  = 

Substituting  these  into  the  Integral  equations,  one  obtains  for  the  velocity  boundary  thick¬ 
ness. 


For  z  =  0 


8  “  « 


/  y  4<li 


/■*  */ 

/  y  dz 


■JT 

20 


5* 

0 


y: 


»  ■  8.0 


-  (f 

The  thermal  boundary  layer  is  given  by 


(125) 

(126) 

(127) 

(128) 


.  izgoIi-sNe,  (-i^)'] 


p 


(129) 


12  ao 


(130) 
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At  the  origin 


For  air  6*  ■  1.92  . 

For  larger  values  of  z 


For  z  ® 


3  0  o 

8«. 


a!  « 


5  N 


5Npr" 


Pr 


8»  =  1.71  Np‘‘ 


12  a  0  V/9Pg 


;;y'4dz 


P 


2  d  z 


8 


T 

■p* 

8i 

8* 

The  heat  transfer  coefficient  is  given  by 

dJ 


h, - 5_ 

8o  '  dy  K-.o 


M  2  Pc  ) 


'/» 


'/z  r  /n  y  ^  d*  1^2 


•  '2  r  7  1 

8«(l2a^)  [ - - — j 


h  s 

X 


k  (  2|»  c,  i 


I  y  J 


(131) 

(132) 

(133) 

(134) 


(135) 


.  Np/' 

8«. 

=  Np," 

(136) 

_ 

=  Np^  ® 

(137) 

(138) 


(139) 
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The  values  for  ( — -■ 

z  >  I 


are  given  by  the  following  approximation: 


/  3  _ 

.936973  [ 

1 

yS  1  * 

(1+z*) 

+  0.844934 

—^2  + 
(1  +  2^^ 

0. 700281 

-  0.482326 

2 

0.  1978  86  ' 
( 

(1+2^)® 

—  0.  029514 

2 

0.  027777 

d+z")* 

in”'  z  ] 


- ^-3  -  0.  542473  - - 

(l+z^)’  (H-z®)^ 


*  +  0.18  2664  - - 


.2.« 


(l  +  z*)^ 


b:.*  ---  } 


(I  +zb 


(140) 


Values  between  z  =  0  and  z  =  1  can  be  obtained  from  the  curve  in  Figure  49.  The  curve 
in  this  region  was  obtained  by  graphical  integration.  At  the  origin  where 

y®  '  20 


(/>  C|  I  )'* 


For  the  cubic  solution  the  following  curves  were  assumed: 


•  (4)' 


.92 


Cubic  Solution 


(141) 


(142) 

(143) 


with  e  =  T  -  T^  as  before.  These  equations  meet  the  boundary  conditions  of: 


y  =  0  ,  u  =  0  , 
y  =  8  ,  u  =  U3 , 


d*u 


dy* 


0 ,  9  =  0, 


dy  ® 


=  0 


du 

dy 


0 


y 


'T  * 


e  = 


±9 

dy 


0 


Substituting  these  into  the  integral  equations,  one  obtains  for  the  velocity  boundary  layer 
thickness 
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140  /  2  ^ 

(tv)  — 


*  '®%6 


dz 


•  >i 


713 


z  —  0 


I  r“’/“ 


•I. 


/is 


Jl-  ^ 

2  44 


Therefore: 


The  thermal  boundary  layer  Is  given  by 


6S*  -S*’ 


540  0 a  Ng,.  (  27g;*IiS^)  (tTT*)  1  -1 

=  \  .  'A  J 


3  ,  2  p.  .72 

(27S*  -  28*  )  (—^}  Y 


/  [  '  +  -f-  t 


sS* - 8* 


27S*  -2S*‘ 


j)  (  ,  +  ,e)  1 


d  z 


540  aa 


T  (278* -  2  8*’) 
For  z  =  0, 


j\y, 


2  dz 


135  aa 


-  i—t 

*)  ^  Pc  > 


^  (27  8*  -2  8** 

The  ratio  of  the  boundary  layer  thicknesses  at  the  origin  is  given  by 


8t 


135  a  a 


P  \'/2 


V  Pc  ' 


2  7  8*  -  28** 


/  6 1  N  (  Pc  > 

X  ( - )  - ^ 

W2  '  o/i 


135 


27  8*  -  28* 


61 

’  72 


-I 


Pr 


16  8*  -216  8*  =  -915  Np/' 
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(147) 

(148) 

(149) 

(150) 


(151) 

(152) 
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For  air,  a  graphical  solution  yields  2.088  for  z  =  0.  As  z  -  for  a  Prandtl  number  of  one, 
6*  approaches  1.0.  The  equation  is 


2781  —  2  S  *  =  25.  I  Npr‘' 


(153) 


The  heat  transfer  coefficient  is  given  by 

k  /  dr 


©0  ^  %  =  o  •  "  2  8^ 


h  -  A.  -A.  (zsoafi  )  ‘  > 

*  ■  2  8 


^4 


r  dz 

®'/i3  ■* 


(154) 


3  I,  \  280  o/t 

2  S^^ 


L  «/,3  J 


(155) 


The  function  is  shown  in  Figure  49.  The  values  in  the  region  z  =  0  to 

z  =  1,  are  obtained  by  graphical  integration.  For  values  of  z  greater  than  one,  the  following 
expression  can  be  used. 


/o  r  ^  dz 
*7.3 


{z  -  0.  9I4I6I  f - —  +  tan^  z  \ 

«Vi3  ^  l-l-z^  ! 


+  0.776  297  — 0.630396 

(l  +  z  )*  (I  +z^)® 


-  -  0.  426638  * 


(l+T*)^ 


-  0.37928  3  ^ — 1  +  0.123729  ^ r+ 0.1 1  421 1  - 


(l+z*)* 


(l  +  z^)* 


(  l  +  z* 


—  0.011352  —  _ -  0.010684 - + 


(l+z*)" 


(l  +  z^)® 


- }  (156) 


At  the  origin 


61  y/z  ip  c,  1)''^^ 


2  2.088  ^  72  *  (a/i  )'/* 

Combined  Thermal  Body  Force  and  Corona  Wind 


(157) 


The  case  of  the  vertical  plate  subjected  to  the  combined  electric  field  and  thermal 
body  forces  was  considered  from  the  integral  equation  approach  also.  The  area  where 
the  corona  wind  velocity  is  approximately  constant  was  considered  in  detail.  The 
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pressure  gradient  term 


_P_ 

Pc 


I 

1+2* 


] 


was  not  included  in  the  initial  analysis.  It 


was  decided  that  the  form  of  the  equations  should  be  Investigated  first  in  the  region  where 
the  pressure  was  approximately  constant.  Based  on  the  outcome  of  this  initial  analysis, 
the  pressure  term  could  be  reconsidered. 

It  was  assumed  that  the  velocity  and  thermal  boundary  layer  thicknesses  were  the  same. 
The  following  velocity  and  temnerature  profiles  were  assumed. 


where 

u^  =  velocity  due  to  the  corona  wind 


(158) 

(159) 


u^  =  velocity  due  to  the  free  convection . 


These  profiles  were  substituted  into  integral  equations.  The  resulting  equations  obtained 
are; 


Momentum: 


I  dS  (  51.  62S  Up*  +9.5  Up  u,  ■«- )  qq  dS  (  l5uo+  2  u,  ) 
105  «IK  ”  24  dx 


Energy: 


r  Tw  “  To  \ 

1  "  / 

3  ' 

^  To  J 

'  8  ' 

8-^8  Ouq+u^j]  =60  a 


(160) 


(161) 


An  Intricate  coupling  exists  and  no  simple  analytic  solution  was  found  for  this  set  of 
equations.  It  would  be  possible  to  determine  a  solution  for  specific  cases  through  numeri¬ 
cal  methods,  but  because  of  the  limitations  of  time  no  further  work  was  expended  on  this 
approach.  Because  of  the  difficulty  experienced  even  with  a  constant  pressure  distribution. 


the  variation  of  pressure 


p 


2 

- ! - y-1  was  not  introduced  into  the  analysis. 

I  +  I  * 


Summary  of  Analysis 


The  analysis  of  the  horizontal  plate  based  upon  the  Von  Karman  integral  approach  appears 
to  have  been  fruitful.  The  magnitude  of  the  heat  transfer  coefficient  is  predicted  to  vary  as 
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ii  and  explicit  relationships  are  available  for  the  variation  of  hx .  Because  of  the  need  for 
assuming  velocity  and  temperature  profiles  in  the  integral  method,  it  cannot  be  expected 
that  the  absolute  magnitude  of  the  heat  transfer  coefficient  will  be  precise.  The  trends, 
however,  should  be  representative  of  the  actual  variations  of  hx  if  the  physical  model 
used  is  a  reasonable  approximation  to  the  actual  test  model.  The  validity  of  the  assump¬ 
tion  of  very  low  space  charge  distortion  of  the  field  equations  can  only  be  ascertained 
after  comparison  with  the  test  data. 

The  analysis  of  the  horizontal  plate  by  means  of  the  thin-jet  approach  appears  to  provide 
an  approximate  technique  for  estimating  the  effects  of  the  corona  discharge  on  heat  trans¬ 
fer.  The  difficulty  of  estimating  the  proper  jet  width,  however,  limits  the  utility  of  this 
method. 

In  the  next  section  a  comparison  will  be  made  between  the  analytical  predictions  and  the 
actual  test  data  for  the  horizontal  plate  case. 


CORRELATION  OF  TEST  DATA  WITH  THE 
ANALYTICAL  SOLUTION 

This  section  will  be  devoted  to  the  comparison  of  the  analytical  solutions  and  the  test 
data  obtained  from  the  interferometer.  The  comparison  is  limited  to  the  case  of  the 
horizontal  plate  with  the  bottom  side  of  the  plate  heated.  The  configuration  considered  is 
the  0.004-inch  wire  located  2  cm  and  5  cm  below  the  plate.  Both  the  integral-equation 
solutions  and  the  two-dim'ensional  thin-jet  solution  will  be  compared. 

In  order  to  compare  the  theory  with  the  test  data,  the  constants  are  evaluated  at  the 
mean  temperature  of  the  air  in  front  of  the  plate.  The  values  of  the  thermal  conductivity, 
density,  and  viscosity  of  air  at  various  temperatures  were  taken  from  Jakob  and  Hawkins 
(Ref.  17). 

The  constant  in  the  equation  for  pressure  pc  =  Ci  i  is  obtained  from  Figure  27,  in 
which  the  data  for  the  0.004-inch  wire  is  plotted  for  various  wire  spacings. 

For 


a  =  2  cm 

c,  =  3.  54  X  10'^  _PiliL 

'  Ml  n 


0  .=  5  cm 

-4  P,  S.f. 

c,  =  3.  80  X  10  - 

'  H-° 

These  values  of  Ci  are  derived  from  the  test  data  on  the  pressure  plate  with  the  wire 
located  vertically  and  the  vertical  dimension  of  the  plate  is  7  inches.  However,  the 
corresponding  dimension  of  the  heated  plate  is  10  inches.  The  current  flowing  to  the 
plate  for  aiy  given  pressure  will  be  proportional  to  the  length  of  the  wire  and  the  plate. 
A  test  run  was  made  to  verify  this.  The  pressure  plate  was  masked  off  with  electrical 
insulating  tape  in  a  series  of  steps  and  pressure-current  measurements  were  taken. 
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The  test  results  shown  in  Figure  50  indicate  that  the  linear  relation  holds.  Correcting  Ci 
for  the  difference  in  plate  size; 

a  =  2  cm 

_4  P.  i.f. 

c,  s  2.48  X  10  ^ - 

fta 

a  =  5  c  m 

-4  P»f. 
c,  =  2.66  X  10  - 

fj.a 

For  a  typical  case  ®  T^eon  =  79“  F, 

V  *  1.686  XIO"^  »q.ft  /sec 
k  =  0.  0«5I  BTU  /  hr  -  ft  -  "F' 
p  =  0.  0736  lbs  /  cu  ft  . 

INTEGRAL  EQUATION  SOLUTION 

The  above  values  are  now  substituted  into  the  equations  for  heat  transfer  coefficient, 
as  follows: 


1.  Equation(139)for  the  linear  profile,  assuming  6*  variation  to  be  small,  becomes; 
For  a  =  2  cm 


For  a  =  5  cm 


(162) 


(163) 


2.  Equation (123) for  the  parabolic  profile,  assuming  6*  variation  to  be  small,  becomes: 
For  a  =  2  cm 

I/.  /  I  r*  dt  "T 

h*  =  0.553  (164) 

For  a  =  5  cm 


=  0. 356i’'‘*  ( 


dz 


(165) 
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3.  Equation (155)for  the  cubic  profile,  assuming  6*  variation  to  be  small,  becomes: 
For  a  =  2  cm 


For  a  =  5  cm 


hx  =  0.479  (— 


r  y  dz  X  2 


'•/IS 


)■ 


TEST  DATA  CORRELATION 


(166) 


(167) 


Typical  solutions  of  heat  transfer  coefficient  for  the  a  =  2  cm  case  are  shown  in  Fig¬ 
ures  51,  52,  and  53.  The  test  data  for  the  corresponding  cases  are  shown  in  Figures  54 
and  55. 

From  a  comparison  of  the  test  data  and  the  theoretical  solutions  at  a  =  2  cm,  it  can  be 
seen  that  the  trends  of  the  solutions  are  good.  The  magnitude  of  h*  differs  depending 
upon  the  profile  considered.  This  had  been  anticipated,  as  was  discussed  in  the  previous 
section.  To  afford  a  more  direct  comparison,  the  theory  and  test  data  are  shown  together 
in  Figure  56,  which  indicates  that  the  trends  are  good.  At  the  end  of  the  plate  the  data  and 
theory  diverge.  A  large  portion  of  this  difference  is  believed  to  be  due  to  the  end  effects 
of  the  plate.  The  comparison  for  a  =  5  cm  is  shown  in  Figure  57.  Both  Figures  56  and  57 
show  that  the  theory  predicts  a  somewhat  more  rapid  decrease  of  hx  as  x  increases  than 
is  shown  by  the  test  data. 

Since  it  appeared  that  the  use  of  the  fixed  value  of  6*  might  introduce  sizeable  variations 
in  hx  as  z  became  large,  a  study  was  undertaken  to  determine  the  explicit  variation  of  6* 
for  all  three  velocity  distributions.  The  variation  of  6*  with  x  is  given  by 


For  the  parabolic  case: 


2 

5*  = 


4  N 


Pr 


8-»-  ( 5  -  S  ■)«■ ) 


For  the  linear  case: 


i  ^  r 

v*''2 


2 

i  *  = 


Npr”' 


f*  1  dz 
!lr '^dz  ^ 


For  the  cubic  case: 


25.05 


- 3  "Pr 

278*  -28* 


yVz  dz 

r 


} 


61 


/»3 


(168) 


(169) 


(170) 
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The  values  of  the  integrals  were  determined  by  graphical  integration  of  the  numerator. 

For  specific  values  of  z,  the  values  of  6*  could  be  determined.  The  linear  case  provided 
an  explicit  solution 

-./ 

x* 

which  exliibits  the  familiar  Np^-1/3  form  found  in  the  forced-convection  theory.  The  values 

of  6*  for  the  parabolic  and  cubic  cases  were  obtained  by  a  careful  graphical  solution  of  the 
4th  order  and  5th  order  equations.  The  variation  of  6*  is  shown  in  Figure  58. 

The  heat  transfer  coefficients  modified  by  the  inclusion  of  the  variation  of  6*  are  shown 
in  Figures  59,  60,  61,  62,  63,  and  64.  The  test  data  for  a  =  5  cm  is  shown  in  Figure  65.  A 
better  comparison  of  the  modified  theory  and  the  test  data  can  be  found  in  Figure  66  where 
the  a  =  2  cm  curves  and  test  data  are  compared  at  two  values  of  the  current.  The  agree¬ 
ment  in  shape  of  the  theoretical  curves  is  now  considered  to  be  excellent.  The  end  effeas 
of  the  plate  are  clearly  evident  in  the  test  data.  Since  the  changes  made  in  6*  had  no  effect 
at  the  origin,  the  small  deviation  between  test  and  theory  remain.  This  difference  is  re¬ 
flected  in  the  relatively  uniform  displacement  of  the  theoretical  curve  from  the  test  data. 
Since  one  of  the  most  sensitive  parameters  to  determine  the  entire  investigation  was  the 
absolute  magnitude  of  pressure,  it  is  not  altogether  unexpected  that  differences  in  magni¬ 
tude  exist.  Figure  67  presents  the  comparison  for  a  =  5  cm.  Similar  results  to  those  pre¬ 
sented  for  the  2  cm  case  are  evident.  Tte  shape  of  the  curves  is  in  good  agreement  except 
where  the  end  effeas  of  the  plate  are  evident.  The  absolute  magnitudes  of  the  heat  transfer 
coefficients  also  differ  somewhat  from  the  test  data. 

The  wide  variation  of  heat  transfer  coefficient  with  corona  wind  over  the  plate  implies 
a  sizeable  variation  of  the  temperature  profile  in  the  boundary  layer.  This  variation  is 
presented  for  a  typical  case  in  Figure  68.  The  temperature  profiles  are  plotted  at  several 
locations  along  the  plate.  The  profiles  do  vary,  and  the  slopes  of  the  curves  increase 
greatly  near  the  origin  under  the  impinging  jet  of  air. 

The  theory  predicted  that  the  heat  transfer  coefficient  should  vary  as  ii  under  the  aaion 
o^  the  corona  wind.  To  verify  that  such  a  trend  exists,  the  values  of  h*  were  plotted  against 
i*.'  The  results  are  shown  in  Figures  69  and  70.  The  current  varies  from  0.32  /la  to  180  u  a 
for  the  2  cm  case,  and  from  2.6  ya  to  180  ua  for  the  5  cm  case.  Data  are  plotted  for  several 
locations  along  the  horizontal  plate.  The  data  at  x  =  0  for  the  a  =  2  cm  case  were  very  er¬ 
ratic  at  high  currents  and  are  not  included.  It  can  be  seen  that  straight  lines  passed  through 
the  points  correlate  with  the  data  quite  well.  Consequently  the  prediction  that  the  heat  trans¬ 
fer  coefficient  should  vary  as  ii  is  substantiated.  Figure  71  presents  data  on  the  variation 
of  the  mean  heat  transfer  coefficient  with  current.  The  ratio  of  the  average  heat  transfer 
coefficient  over  the  surface  of  the  plate  with  applied  field  to  the  average  free  convection 
heat  transfer  coefficient  is  plotted  against  ii.  It  can  be  seen  that  a  linear  relationship 
also  would  hold  well  for  this  case. 

It  is  apparent  from  a  review  of  the  theory  and  the  experimental  work,  that  large  changes 
in  heat  transfer  over  the  free-convection  case  can  be  achieved  through  the  use  of  the  corona 
wind.  Substantial  increases  occur  for  very  small  currents.  The  power  required  for  the 
action  is  low;  some  typical  values  are  shown  in  Table  1.  These  values  are  for  some  of  the 
curves  shown  in  Figures  54,  55,  and  65. 
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Table  1 


Power  Requirements  for  Corona  Wind,  Horizontal  Plate.  0.004-Inch 
Wire  at  Two  Spacings. 


Wire  Spacing 
(cm.) 

Current 

Voltage 

(KV) 

Power 
(watts  ) 

2 

0 

0 

0 

2 

.32 

5.00 

.0016 

2 

1.05 

5.00 

.0053 

2 

2.54 

5.17 

.0131 

2 

4.55 

5.67 

.0258 

2 

0 

0 

0 

2 

5.6 

6.0 

.0336 

2 

10.2 

6.34 

.0648 

2 

19.1 

6.67 

.127 

2 

41.0 

8.00 

.328 

2 

79.5 

8.83 

.703 

5 

0 

0 

0 

5 

2.6 

8.33 

.022 

5 

4.55 

8.67 

.039 

5 

10.2 

10.00 

.102 

5 

21.2 

11.85 

.251 

5 

39.5 

14.19 

.561 

5 

78.5 

18.18 

1.429 

- 

These  data  correspond  to  the  data  in  Figures  54.  55,  and  65. 
THIN  JET  SOLUTION 


T^e  last  correlation  to  be  made  is  between  the  thin  jet  and  the  test  data.  The  equations 
ied  were; 

Stagnation  region  (Ref.  10): 

^  4u 

X  <  -  b 

ttU 

h,  =  0.570  k 


I  \  2 

•  ~  )  (172) 
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h 


X 


k 


(0.  332  Np 


(173) 


Study  of  the  pressure  distribution  revealed  that  an  approximate  equivalent  jet  half-width 
for  a  *  2  cm  would  be  b  -  a/2.  Using  this  value  and  substituting: 


Pj.  =  c,  I  2  2.478  X  10*^  i  (174) 


The  equation  for  the  stagnation  region  becomes: 


hx 


(0.5  70  Np“  '*  )  k 


I  2  p  v'/a  /  2  X  2.4781  ^'*^4 
\o  fj.  )  \  p  ) 


Typically  for  air  at  79*F  and  a  =  2  cm: 

'/ 

=  I.  98i  ♦ 


(175) 


(176) 


This  expression  holds  from  x  =  0,  to  x  =  2a/  .  Beyond  x  =  2a/  the  Blasius  solution  is 
used.  Fitting  the  Blasius  solution  so  that  at  x  =  2a/  the  heat  transfer  coefficients  of 
the  two  regions  are  the  same,  results  in  the  following  expression  when  x  >  2a/w: 


where  x  is  given  in  centimeters. 


2.23  l 
x'/* 


(177) 


These  relationships  are  plotted  in  Figure  72  for  a  =  2  cm.  Comparison  of  this  curve 
with  the  test  data  of  Figure  55  indicates  a  fair  degree  of  correlation.  The  right  order 
of  magnitude  is  achieved  and  the  shapes  of  the  curves  approximate  those  for  the  test 
data  except,  of  course,  in  the  region  of  the  ends  of  the  plate.  The  primary  problem  in 
this  method  lies  in  the  choice  of  a  suitable  jet  half-width,  b.  No  accurate  nor  logical 
method  presents  itself  by  which  to  determine  b.  Since  both  magnitude  and,  to  some  ex¬ 
tent,  shape  are  affected  by  the  choice  of  b,  this  method  can  be  considered  only  an  approx¬ 
imation. 


EXPLORATORY  TESTS  ON  FREE  CONVECTION 
WITH  OTHER  ELECTRODES 

In  the  over-all  investigation  of  the  electrostatic  effects  on  free-convection  heat  trans¬ 
fer,  many  electrode  configurations  were  used.  At  the  outset  of  the  testing,  a  Model  “T” 
Ford  automotive  spark  coil  was  used  as  a  high-voltage  source.  The  voltage  wave  form 
was  found  to  be  very  erratic,  including  both  a  D.C.  component  and  random  A.C.  com¬ 
ponents.  The  spark  coil  was  useful  in  providing  some  interesting  pictures  of  effects  of 
the  corona  discharge  on  the  interferometer  patterns.  Results  of  tests  made  with  a  two- 
dimensional  row  of  pins  and  the  spark  coil  are  shown  in  Figure  73.  View  A  illustrates 
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the  effect  on  the  isothermal  lines  an  instant  after  the  field  is  applied;  a  large  swirl  can 
be  seen  moving  upward  at  the  top  of  the  picture.  View  B  shows  the  effect  of  the  corona 
discharge  under  continuous  application  of  the  field.  In  both  pictures  the  row  of  needles 
is  represented  by  the  gray  line  pointing  upward  towards  the  center  of  the  plate.  The  dual 
lines  on  the  left  side  of  the  picture  are  the  support  rods  for  the  needles. 

Figure  74  shows  the  results  of  an  initial  test  in  which  a  grid  of  0.010-inch-diameter 
wires  was  placed  in  front  of  the  heated  plate.  (The  wire  grid  is  shown  in  the  rear  of  Fig¬ 
ure  7.)  The  electric  field  is  impressed  between  the  fine  wires  and  the  heated  plate.  When 
voltagp  is  applied  by  means  of  the  automotive  spark  coil,  the  thermal  boundary  layer  in 
the  vicinity  of  the  wires  is  pulled  out  decidedly.  These  pictures  were  taken  in  the  initial 
exploratory  phase  of  the  work.  The  photography  is  very  poor,  but  the  pictures  are  in¬ 
cluded  because  the  phenomenon  they  show  had  an  important  impact  on  subsequent  testing. 
Many  of  the  subsequent  vertical  plate  tests  with  a  single  wire  described  in  a  previous 
section  were  aimed  at  understanding  this  pullout.  The  pullout  phenomenon  appears  to  in¬ 
dicate  an  action  counter  to  that  of  corona  wind,  and  even  counter  to  the  action  which  would 
be  produced  by  dielectric  streaming.  As  indicated  in  Reference  1,  dielectric  streaming 
would  occur  as  the  result  of  a  body  force  produced  by  the  nonuniform  electric  field.  How¬ 
ever,  the  density  gradient  and  thus  the  dipole  moment  per  unit  volumejis  in  the  direction 
to  cause  the  dielectric  streaming  to  move  away  from  the  fine  wires.  As  discussed  under 
testing,  it  is  believed  that  the  pullout  occurs  as  an  interaction  between  fluid  streams,  and 
one  of  these  streams,  the  boundary  flow,  is  highly  viscous.  Thus,  it  appears  that  the  first 
impression  of  the  picture  may  be  misleading. 

Figure  75  illustrates  the  results  of  tests  with  the  same  wire  grid  but  using  the  high 
voltage  provided  by  the  D.C.  power  supply.  At  the  very  low  currents,  between  1  oa  and 
5  ua,  the  same  pullout  phenomenon  occurred  as  was  observed  with  the  spark  coil;  pictures 
in  the  low  current  range  are  not  included.  View  A  is  a  reference  with  no  field  applied.  The 
remaining  pictures  indicate  the  effect  of  increased  currents  on  the  isothermal  lines.  The 
thermal  boundary  layer  is  made  thinner  by  the  application  of  the  electric  field,  and  the 
entire  flow  field  becomes  erratic  and  apparently  very  turbulent. 

Test  results  using  a  specially  prepared  surface  for  the  heated  plate  are  shown  in  Fig¬ 
ure  76.  This  surface  consisted  of  a  series  of  0.(X)4-inch  wires  attached  to  the  surface  of 
an  insulating  sheet  of  plastic.  The  wires  were  spaced  one  inch  apart^parallel  to  each  other 
and  horizontal  to  the  plastic  sheet,  which  was  cemented  to  the  front  of  the  heated  plate. 
The  fine  wires  were  located  away  from  the  plate  surface  facing  the  free  room  air.  Drop¬ 
ping  resistors  were  placed  between  each  of  the  wires,  and  were  so  chosen  that  approxi¬ 
mately  a  current  of  10  pa  per  wire  would  flow  for  an  applied  voltage  of  30  KV.  It  was 
believed  that  this  configuration  would  provide  for  local  sources  of  ions  within  the  bound¬ 
ary  layer.  By  placing  the  other  electrode  at  various  locations  and  applying  a  voltage,  the 
possible  effects  on  the  boundary  layers  could  be  observed  with  the  interferometer.  View 
A  illustrates  the  reference  case  with  no  applied  field.  The  wires  are  so  small,  they  cannot 
be  seen  at  the  surface  of  the  plate.  The  external  electrode  is  positioned  1-5/8  inches  be¬ 
low  the  bottom  edge  of  the  plate.  This  electrode,  the  oval  rod  shown  in  the  front  of  Fig¬ 
ure  7,  is  not  visible  in  the  interferometer  pictures.  View  B  shows  the  effects  when  the 
field  is  applied.  Even  at  high  currents  very  little  effect  of  the  field  is  found.  It  is  appar¬ 
ent  that  any  attempt  to  control  fluid  behavior  at  the  surface  of  a  body  through  the  applica¬ 
tion  of  electric  fields  will  have  to  be  done  with  a  great  deal  of  careful  forethought. 

Figure  77  illustrates  the  effect  of  a  uniform  electric  field  on  the  thermal  boundary  layer 
of  the  vertical  heated  plate.  A  flat  6-inch-by- 10-inch  aluminum  electrode,  unheated,  was 
placed  1.1  cm  in  front  of  the  heated  plate.  Views  A  and  B  show  results  before  and  after  a 
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strong  field  was  applied.  No  corona  discharge  occurred/io  measurable  currents  were 
observed,  and  no  changes  to  the  fringe  patterns  are  apparent.  Tests  of  this  type  were 
repeated  using  several  electrode  spacings  and  applied  voltages  up  to  and  including  that 
necessary  for  breakdown  (sparking).  No  effects  on  the  isothermal  lines  were  visible  for 
applied  fields  below  breakdown  potentials.  It  thus  appears  that  application  of  a  uniform 
field  in  air  with  very  low  ionization  produces  no  measurable  changes  on  the  free-convection 
heat  transfer. 

EXPLORATORY  TESTS  WITH  DIELECTROPHORESIS 

This  section  will  cover  very  briefly  a  series  of  tests  made  to  explore  the  body  forces 
obtainable  through  the  use  of  the  dielectrophoresis  body  force.  The  phenomenon  of  dielec¬ 
trophoresis  is  explained  in  detail  in  Reference  1.  In  summary,  it  is  a  force  which  can  act 
on  a  particle  placed  in  a  nonuniform  electric  field  even  if  that  particle  does  not  possess 
a  net  electrical  charge.  By  virtue  of  the  charge  displacements  occurring  on  bodies  in 
electric  fields,  an  electrical  body  force  arises.  Experiments  conducted  by  Gemant  and 
Pohl  (Refs.  22.  23)  have  indicated  that  motion  of  the  particles,  either  solid  or  liquid,  can 
be  achieved.  To  verify  these  results,  similar  tests  were  conducted.  In  addition  some  orig¬ 
inal  work  was  daie  on  the  motion  of  drqslets  of  water  in  a  simulated  “Zero-G”  model. 

DESCRIPTION  OF  TEST  APPARATUS 

Figure  78  shows  the  entire  test  set-up.  At  the  right  is  the  D.C.  high-voltage  power  sup¬ 
ply;  next  to  it  is  a  neon  transformer  ctxmected  to  a  Variac  to  provide  a  variable  source  of 
high  A.C.  voltage.  The  test  model,  located  at  the  center  of  the  picture,  is  shown  enlarged 
in  Figure  79.  Several  electrodes  of  various  shapes  are  shown  at  the  left  of  the  glass  con¬ 
tainer.  The  test  model  is  a  glass  container  approximately  8  inches  in  diameter;  one  elec¬ 
trode  lead  extends  through  a  sealed  stem  at  its  right  side.  For  these  tests,  one  of  the 
electrodes  was  fastened  to  the  lead  in  the  container,  and  another  was  mounted  in  the  elec¬ 
trode  holder  directly  above  the  center  ctf  the  container.  The  electrical  field  was  applied 
between  the  upper  and  the  lower  electrodes.  The  container  was  filled  with  fluid  consisting 
of  a  mixture  of  one  half  each  of  CCl*  and  silicone  D.C.  200  fluid.  The  fluids  were  mixed 
in  an  attempt  to  obtain  a  fluid  that  was  nonconducting  but  had  the  same  specific  gravity 
as  that  of  water.  It  had  been  anticipated  that  the  fluids  were  completely  miscible,  but  an 
interface  formed  between  two  phases  of  the  fluid  at  which  it  was  possible  to  suspend  drop¬ 
lets  of  water. 

TEST  RESULTS 

Figures  80  and  81  illustrate  the  electric  wind  effects.  The  point  used  as  the  central 
electrode  is  located  just  above  the  surface  of  the  liquid.  The  lower  electrode,  directly 
beneath  it,  is  in  the  fluid  at  the  bottom  of  the  container.  When  a  field  is  applied,  waves 
are  generated  on  the  surface  of  the  liquid  as  the  corona  wind  strikes  the  surface.  This 
phenomenon  is  the  same  as  that  considered  previously  in  the  sections  on  heat  transfer. 

The  next  series  of  pictures  indicate  the  effect  of  a  nonuniform  A.C.  field  on  the  fluid. 

The  central  elearode,  shown  in  Figure  82,  is  a  round  ball  approximately  1/8  inch  in  dia¬ 
meter  located  3/16  inch  above  the  surface  of  the  fluid.  The  lower  electrode  is  the  same 
as  was  used  before.  When  the  A.C.  voltage  is  applied,  the  fluid  jumps  up  and  clings  to  the 
surface  of  the  ball,  as  shown  in  Figure  83.  A  ring  central  electrode  was  placed  above  the 
surface  as  shown  in  Figure  84.  Figure  85  shows  the  fluid  clinging  to  the  wire  under  the 
action  of  the  field.  These  experiments  merely  represent  pictorial  verification  of  Gemant’s 
and  Pohl’s  work. 
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The  last  pictures  illustrate  the  action  of  the  nonuniform  field  on  suspended  water  drop¬ 
lets.  Figure  86  shows  the  fluid  model.  The  upper  electrode  is  a  small  round  surface  elec¬ 
trically  insulated  by  a  very  heavy  coating  of  wax,  and  the  lower  electrode  is  a  sheet  of  alu¬ 
minum  foil  placed  below  and  outside  the  glass  container.  The  nonuniform  field  is  established 
between  these  two  electrodes.  Two  relatively  large  droplets  of  water  can  be  seen  suspend¬ 
ed  upon  the  interface  existing  within  the  fluid  model.  After  the  field  was  applied,  the  water 
droplets  moved  slowly  toward  the  central  electrode  and,  upon  reaching  it,  coalesced  to 
form  a  single  droplet.  This  sequence  portrays  the  actual  existence  of  the  nonuniform  field 
body  force.  Considerable  care  had  to  be  exercised  in  conducting  the  tests  in  order  to  ob¬ 
serve  the  phenomenon.  Whenever  a  bare  central  electrode  was  used,  the  droplet  would  be 
repelled  and  broken  up  upon  reaching  it. 

These  exploratory  tests  are  Included  solely  to  demonstrate  that  nonuniform  field  effects 
actually  exist.  It  was  not  possible  in  the  time  available  to  obtain  and  use  accurate  instru¬ 
mentation.  Approximately  4000  volts  A.C.  were  required  to  move  the  droplets  of  water 
but,  unfortunately,  currents  were  not  measurable  with  the  A.C.  milliammeter  available. 

The  pictures,  however,  do  indicate  that  some  measure  of  fluid  orientation  and  control 
are  possible  with  the  nonuniform  electric  fields. 


SUMMARY  OF  ELECTRIC  FIELD-FLUID  INVESTIGATION 

Starting  with  the  fundamental  question,  what  are  the  possible  effects  of  an  electrostatic 
field  on  the  boundary  layer  and  heat  transfer  in  a  fluid,  an  entire  field  of  investigation  was 
opened.  The  review  of  literature  and  theory  in  several  fields  of  endeavor  revealed  that  a 
broad  range  of  influences  were  possible.  The  question  that  followed  was,  were  the  phe¬ 
nomena  controllable  and  useful? 

Two  series  of  experiments  were  then  conducted  to  determine  whether  or  not  the  influences 
could  be  controlled  or  made  to  provide  useful  results.  One  series  consisted  of  determining 
the  effects  of  corona  wind  on  free -convection  heat  transfer.  An  extensive  experimental  and 
analytical  program  proved  that  large  changes  in  heat  transfer  occurred  with  the  application 
of  electric  fields.  The  correlation  between  theory  and  test  data  appeared  to  be  very  good. 
While  the  actual  magnitudes  of  the  heat  transfer  coefficients  were  not  predicted  exactly, 
the  indicated  trends  were  in  excellent  agreement.  The  second  series  of  experiments  were 
purely  exploratory  and  were  accomplished  to  demonstrate  some  actions  of  the  dielectro¬ 
phoresis  phenomenon. 


CONCLUSIONS  AND  RECOMMENDATIONS 

As  a  result  of  the  over-all  survey  to  determine  the  possible  interactions  between  elec¬ 
tric  fields  and  fluids,  it  was  concluded  that  many  possible  actions  exist.  The  actions 
include: 

1.  Ion-drift  phenomena  due  to  an  applied  field 

2.  Charged-particle  motion  in  an  A.C.  field 

3.  Phoresis  motions 

4.  Dielectrophoresis  phenomena 

5.  Changes  in  fluid  properties  with  applied  fields. 
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Because  of  these  actions,  body  forces  are  available  of  sufficient  magnitude  to  perform 
useful  functions.  Prior  experiments  have  indicated  that  both  partially  ionized  and  neutral 
fluids  can  be  influenced.  It  appears  logical  to  utilize  these  body  forces  in  carefully  select¬ 
ed  regions  of  the  fluid,  in  the  boundary  layer,  for  example,  to  provide  changes  in  heat  trans¬ 
fer  or  drag. 

The  experimental  work  conducted  in  this  study  clearly  indicates  that  a  sizeable  change  in 
heat  transfer  can  be  obtained  through  the  use  of  the  corona  discharge.  The  application  of 
the  Von  Karman  integral  equations  and  electric  field  theory  lead  to  an  analytical  solution 
which  correlates  well  with  the  test  data  for  a  typical  test  configuration.  The  demonstration 
tests  of  dielectrophoresis,  likewise,  indicate  that  motion  of  neutral  fluids  and  particles 
is  possible  with  nonuniform  A.  C.  electric  fields. 

It  is  believed  that  the  phenomena  involved  in  electrofluldmechanics  can  have  significant 
importance  in  a  wide  range  of  fluid  mechanic  and  heat  transfer  problems.  Much  basic 
work  must  be  done,  however,  to  understand  the  various  phenomena,  their  limitations, 
and  their  scope. 

It  is  recommended  that  a  continuing  effort  be  expended  in  this  area.  Many  possible  ave¬ 
nues  of  research  are  evident,  including  effects  on  viscosity,  thermal  conductivity,  trans¬ 
ition,  separation,  drag,  and  heat  transfer.  Potential  areas  of  application  may  be  large  if 
the  phenomena  can  be  controlled  at  will. 
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Figure  3.  Close-up  of  camera  mounting  on  interferometer 
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Typical  interferometer  installation  of  heated  plate  and  electrode 


ASD-TDR-62-650 


60 


i 


Figure  7.  Electrodes  for  interferometer  tests 
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Figure  8.  Over-all  view  of  inclined  manometer  showing  spotlight  and  7  x  1,7-inch  plate 
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Over-all  view  of  pressure  survey  arrangement  using  Chattock  micromanometer 
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A.  No  electric  field 


Figure  IJ.  Interferometer  pictures  of  the  effects  of  corona  discharge  on  heat  transfer.  Vertical 
heated  plate,  plate  temperature  approximately  119°F,  room  temperature  690F,  v/ire 
0.004  inch  in  diameter,  2  cm  in  front  of  plate  at  centerline 
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5,9  microamps 


22. 8  mior Damps 


O.  10.0  KV  74.2  microamps  U.  U.lKV  195  microamps 

Figure  13  (Coni'  d).  Interferometer  pictures  of  the  effects  of  corona  discharge  onheat  transfer.  Vertical 
heated  plate,  plate  temperature  approximately  II9OF,  room  temperature  6 9°F,  wire 
0.004  inch  in  diameter,  2  cm  in  front  of  plate  at  eenterline 
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7  KV  19.9  microamps  D.  11.7  KV  162  microamps 

Figure  15.  Effects  of  corona  discharge  on  heat  transfer.  Vertical  heated  plate, 
room  temperature  64°F,  AT  held  constant  at  Z7'^F,  0.004-inch  wire 
located  2  cm  from  plate  at  centerline 
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16.  Heat  transfer  coefficients  for  vertical  plate,  0.004- inch  wire  located 
2  cm  from  centerline,  corresponding  to  photographs  of  Figure  15 
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1-5KV  19.9  microamps  D.  23. 3  KV  158  xnicroamps 

Fig-ure  17.  Effects  of  corona  discharge  on  heat  transfer.  Vertical  heated  plate, 
room  temperature  hl^F,  AT  held  constant  at  34°F,  0.004- inch  wire 
located  5  cm  from  plate  at  centerline 
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Figure  19.  Effects  of  corona  discharge  on  heat  transfer.  Vertical  heated  plate. 

room  temperature  53°F,  i^T  held  constant  at  34'’f,  0.004-lnch  wire 
located  10  cm  from  plate  at  centerline 


Figure  21.  Effects  of  corona  discharge  on  heal  transfer.  Horizontal  heated  plate, 
room  temperature  60°F,  held  constant  at  27°F,  0.004-inch  wire 
located  2  cm  from  plate  at  centerline 
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Figure  23.  Effects  of  corona  discharge  on  heat  transfer.  Horizontal  heated  plate 
room  temperature  61®F,  AT  held  constant  at  34oF,  0.004-inch  wire 
located  2  cm  from  plate  at  centerline 
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Figure  42.  End  view  of  horizontal  plate  with  bottom  heated,  subjected  to  corona 
discharge  In  air 


Figure  43.  End  view  of  vertical  plate  with  front  heated,  subjected  to  corona 
discharge  in  air 
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Figure  46.  Comparison  of  the  narrow  ideal  jet  pressure  distribution  with  the  data  of  corona 
discharge  from  a  0.004-inch  wire  with  Z  cm  spacing 
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Figure  51.  Theoretical  variation  of  heat  transfer  coefficient,  parabolic  profile, 
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Figure  52.  Theoretical  variation  of  heat  transfer  coefficient,  linear  profile, 
a  =  2  cm 
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Figure  53.  Theoretical  variation  of  heat  transfer  coefficient,  cubic  profile, 
a  =  2  cm 
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Figure  56.  Direct  comparison  d  theory  and  test  data,  a  =  2  cm, 
wire  diameter  0.004  inch 
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Figure  58.  Variations  of  6#  for  various  profiles 
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Figure  59.  Modified  heat  transfer  coefficients,  parabolic  profile,  a  =  2  cm 
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Figure  66.  Direct  comparison  of  modified  theory  with  test  data  for  a  =  2  ci 
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horizontal  plate  under  Influence  of  corona  discharge,  a  =  2  cm, 
current  5, 6  microamperes 


Figure  69.  Variation  of  local  heat  transfer  coefficient  with  corona  current, 
a  =  2  cm,  wire  diameter  0.004  inch 
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Variation  of  local  heat  transfer  coefficient  with  corona  current, 
a=  5  cm,  wire  diameter  0,004  inch 
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L  of  mean  heat  transfer  coefficient  for  horizontal  plate, 
3h  wire 


A.  Just  after  field  applied  B.  Steady  state  with  field 

Figure  73.  Corona  di.scharge  to  a  heated  plate  using  a  row  oi  needles  and  spark  oo: 


A.  Approximately  1  cm  spacing 


B.  Approximately  2  cm  spacing 


Figure  74.  Electrical  field  applied  to  a  wire  grid  located  in  front  of  the  vertical 
heated  plate  using  a  spark  coil 


10.2  KV  19.8  microamps  D.  12. 5  KV  200  micro  amps 

Figure  75.  Interferometer  pictures  of  the  rffects  of  corona  discharge  on  heat  transfer. 
Vertical  heated  plate,  room  temperature  52°F,  A  T  held  constant  at  340B', 
0.010-inch  wire  grid  located  horizontally  at  2  cm  from  plate,  D.C.  voltage 
applied 
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B.  39.2  KV  Current  greater  than  50  microamps 
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Figure  76.  Interferometer  pictures  of  the  effects  of  corona  discharge  on  heat  transfer 
Boundary  layer  grid  glued  to  vertical  heated  plate,  plate  temperature 
approximately  100°F,  room  temperature  54'’F,  oval  positive  electrode 
1-5/ 8-lnch  below  bottom  edge  of  heated  plate. 


Zero  Voltage  0  microamps  B.  14.7  KV  0  microamps 

Figure  77.  Effeet  of  uniform  electric  field  on  free  convection  heat  transfer 
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Figure  79.  Close  up  ol  dielectrophoresis  equipment  showing  electrodes 
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Figure  82.  Dielectric  attraction  test,  no  applied  field,  sphere  located  just  above 
liquid  surface 


Figure  83.  Dielectric  attraction  test  with  applied  field 
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Figure  84.  Dielectric  attraction  test,  no  applied  field,  ring  electrode  just 
above  surface 


Figure  85.  Dielectric  attraction  test  with  applied  field 


Figure  86.  Dielectric  influence  on  water  droplets  suspended  in  a  nonconducting 
liquid,  no  applied  field 


Figure  87.  Dielectric  influence  on  water  droplets  suspended  in  a  nonconducting 
liquid,  alter  the  application  of  the  electric  field 
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APPENDIX  I 

CALCULATION  OF  HEAT  TRANSFER  COEFFICIENTS 

The  calculation  of  the  heat  transfer  coefficients  from  the  interferograms  is  a  straight¬ 
forward  process  (Ref.  1)*.  It  is  based  upon  the  fundamental  equations  of  heat  transfer 
and  the  unique  characteristics  of  the  interferometer.  Consider  the  basic  heat  transfer 
equations  first.  In  convective  heat  transfer  from  a  wall,  it  is  assumed  that  the  wall  is  at 
rest.  It  then  follows  that  the  heat  must  pass  from  the  wall  through  the  fluid  layer  by  con¬ 
duction,  and  the  simple  equation  of  heat  conduction  can  be  used. 


dy  y  so 

where  A  is  the  area  through  which  the  heat  passes.  The  defining  equation  for  the  heat 
transfer  coefficient  is: 


a  =  h  A  (  -  To  )  . 


Equating  these  two  expressions,  one  finds; 

k  dT  1 

h  =  - - 

-  To  d  y  yr  0 


The  value  of  k  is  determined  at  a  temperature  representative  of  the  fluid  layers  adjacent 
to  the  wall.  Following  Sparrow  (Ref.  2),  the  reference  temperature  used  in  evaluating  the 
fluid  properties  is 

Tr  *  T*  -  0.  38  (  T^  -  To  ) 


The  values  of  the  properties  of  air,  including  density,  viscosity,  and  Prandtl  number  were 

taken  from  References  1  and  3.  -*4^^ — I  is  determined  from  the  interferometer  data. 

ly  =  o 

Eckert  and  Soehngen  have  derived  an  equation  for  the  temperature  difference  between 
ambient  and  any  fringe  lathe  boundary  layer  (Ref.  4).  A  complete  tabulation  of  values 
of  this  equation  is  given  by  Evalenko  (Ref.  1),  and  these  values  are  reproduced  in  Table 
2.  To  use  the  table,  it  is  necessary  to  know  the  density,  and  the  number  of  the  fringe 
being  considered.  The  fringe  temperature  is  then  given  by 


Xq  e 
Lc  yo 


In  order  to  determine  the  temperature  gradient  at  the  wall,  it  is  necessary  to  measure 
the  spacing  of  the  fringes  precisely.  The  interferometer  pictures  were  read  by  means  of 
a  Gaertner  comparator.  This  instrument  can  read  distances  to  within  0.0001  inch.  A  small 
“T”  is  included  at  the  bottom  of  each  photograph  to  provide  a  reference  length  of  three 
inches  from  which  a  scale  factor  could  readily  be  obtained  for  each  photograph.  Lines 
were  carefully  scribed  on  the  photographs  at  selected  stations  along  the  heated  plate. 


*Numbers  in  parentheses  refer  to  references  for  Appendix  I. 
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Table  2 


Values  of  the  Properties  of  Air, 


) 


F 

p  « .0710 

Diff 

p  -  .0720 

Diff 

p  « .0730 

Diff 

p  =  .0740 

1 

.00840 

12 

.00828 

11 

.00817 

12 

.00805 

2 

.01694 

25 

.01669 

23 

.01646 

22 

.01624 

3 

.02564 

38 

.02526 

36 

.02490 

34 

.02456 

4 

.03447 

51 

.03396 

48 

.03348 

47 

.03301 

5 

.04345 

64 

.04281 

60 

.04221 

60 

.04161 

6 

.05260 

79 

.05181 

74 

.05107 

71 

.05036 

7 

.06192 

95 

.06097 

87 

.06010 

85 

.05925 

8 

.07140 

111 

.07029 

100 

.06929 

100 

.06829 

9 

.08104 

126 

.07978 

115 

.07863 

115 

.07748 

10 

.09087 

139 

.08948 

135 

.08813 

128 

.08685 

11 

.10087 

159 

.09928 

147 

.09781 

144 

.09637 

12 

.11106 

178 

.10928 

162 

.10766 

162 

.10604 

13 

.12145 

197 

.11948 

179 

.11769 

178 

.11591 

14 

.13202 

216 

.12986 

195 

.12791 

196 

.12595 

15 

.14278 

233 

.14045 

215 

.13830 

213 

.13617 

16 

.15377 

254 

.15123 

233 

.14890 

233 

.14657 

17 

.16496 

276 

.16220 

252 

.15968 

250 

.15718 

18 

.17639 

298 

.17341 

272 

.17069 

273 

.16796 

19 

.18804 

323 

.18481 

291 

.18190 

292 

.17898 

20 

.19990 

344 

.19646 

315 

.19331 

312 

.19019 

21 

.21201 

367 

.20834 

339 

.20495 

332 

.20163 

22 

.22437 

396 

.22041 

357 

.21684 

358 

.21316 

23 

.23699 

421 

.23278 

382 

.22896 

382 

.22514 

24 

.24987 

448 

.24539 

409 

.24130 

406 

.23724 

25 

.26303 

474 

.25828 

437 

.25391 

431 

.24960 

26 

.27645 

506 

.27139 

461 

.26678 

456 

.26222 

27 

.29018 

539 

.28479 

488 

.27991 

484 

.27507 

28 

.30418 

569 

.29849 

516 

.29333 

513 

.28820 

29 

.31852 

605 

.31247 

547 

.30700 

541 

.30159 

30 

.33317 

640 

.32677 

576 

.32101 

574 

.31527 

31 

.34813 

672 

.34141 

612 

.33529 

605 

.32924 

32 

.36344 

712 

.35632  . 

642 

.34990 

641 

.34349 

33 

.37908 

747 

.37161 

678 

.36483 

676 

.35807 

F  is  the  fringe  number 


p  is  the  density  of  the  air  based  upon  T- 


124 


ASD-TDR-62-650 

Distance  readings  were  made  for  all  fringes  at  a  given  station,  and  corrected  by  applying  the 
scale  factor  for  that  photograph.  Knowing  the  temperature  of  each  fringe  and  its  spacing, 
it  is  then  possible  to  plot  the  temperature  profiles.  From  the  temperature  profiles,  the 
slope  at  the  wall  was  determined  for  several  typical  cases.  For  the  same  cases  the  slopes 
based  upon  the  inner  two  fringes  were  determined  and  compared  with  the  slopes  based 
upon  the  plotted  temperature  profiles.  Since  the  two  slopes  agreed  quite  well,  the  deter¬ 
mination  of  the  remaining  temperature  gradients  was  based  upon  the  slopes  determined  by 
the  temperatures  and  spacings  of  the  inner  two  fringes. 
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APPENDIX  II 

LIST  OF  EQUIPMENT  USED 

The  nomenclature  for  the  equipment  used  in  the  experimental  phase  is  listed  in  this 
section. 

1.  High  voltage  D.C.  power  supply 

Type:  D.  C.  power  supply 
Rating:  0-2  ma 
0-50  KV 

Make:  Sterling  Instruments  Co.,  Detroit,  Michigan 
Model:  PS-50 

2.  Microammeter 

Type:  Illuminated  Spot  microammeter 

Rating:  0-1-3-10-100-300-1000  ua 

Make:  Greibach  Instruments  Corp.,  New  Rochelle,  N.  Y. 

Model:  Number  500 

3.  High  Voltage  Voltmeter 

Type:  Electrostatic  voltmeter,  double  pivoted  moving  vane 
Rating:  0-10-20-30-40  KV 

Make:  Sensitive  Research  Instruments  Corp.,  New  Rochelle,  N.  Y. 

Model:  ESH 

4.  High  Voltage  Voltmeter 

Type:  Vacuum  tube  voltmeter  plus  high  voltage  probe 
Rating:  0-50  KV 

Make:  Radio  Corporation  of  America,  Harrison,  N.  J. 

Model:  WV-97A  and  WG-289  High  Voltage  Probe 

5.  Thermocouple  Potentiometer 

Type:  Portable  precision  potentiometer 

Make:  Rubicon  Instruments  Div.,  Minneapolis-Honeywell  Regulator  Co., 
Minneapolis,  Minn. 
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Inclined  Minometer 


T)rpe:  40  tube  band  with  reservoir 

Make;  Ann  Arbor  Instrument  Works,  Ann  Arbor,  Michigan 
Specific  gravity  of  fluid  .956 

Micrpmanometer 

T5rpe;  ChattocI:  two  fluid  micromanometer  (Kerosene  and  water) 
Make;  Frank  G.  Wahlen,  Chicago,  Ill. 

Traversing  Microscope 

Type;  micrometer  traverse  -  Range  0-4  '  traverse 
Make;  Gaertner  Scientific  Corp.,  Chicago,  Ill. 

Model;  Gaertner  Comparator 


